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Abstract. We review the theory of continuous-variable entanglement with special emphasis 
on foundational aspects, conceptual structures, and mathematical methods. Much attention 
is devoted to the discussion of separability criteria and entanglement properties of Gaussian 
states, for their great practical relevance in applications to quantum optics and quantum 
information, as well as for the very clean framework that they allow for the study of the 
structure of nonlocal correlations. We give a self-contained introduction to phase-space and 
symplectic methods in the study of Gaussian states of infinite-dimensional bosonic systems. 
We review the most important results on the separability and distillability of Gaussian states 
and discuss the main properties of bipartite entanglement. These include the extremal 
entanglement, minimal and maximal, of two-mode mixed Gaussian states, the ordering of two- 
mode Gaussian states according to different measures of entanglement, the unitary (reversible) 
localization, and the scaling of bipartite entanglement in multimode Gaussian states. We then 
discuss recent advances in the understanding of entanglement sharing in multimode Gaussian 
states, including the proof of the monogamy inequality of distributed entanglement for all 
Gaussian states. Multipartite entanglement of Gaussian states is reviewed by discussing 
its qualification by different classes of separability, and the main consequences of the 
monogamy inequality, such as the quantification of genuine tripartite entanglement in three- 
mode Gaussian states, the promiscuous nature of entanglement sharing in symmetric Gaussian 
states, and the possible coexistence of unlimited bipartite and multipartite entanglement. We 
finally review recent advances and discuss possible perspectives on the qualification and 
quantification of entanglement in non Gaussian states, a field of research that is to a large 
extent yet to be explored. 
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1. Prologue. 

About eighty years after their inception, quantum mechanics and quantum theory are still an 
endless source of new and precious knowledge on the physical world and at the same time 
keep evolving in their mathematical structures, conceptual foundations, and intellectual and 
cultural implications. This is one of the reasons why quantum physics is still so specially 
fascinating to all those that approach it for the first time and never ceases to be so for those 
that are professionally involved with it. In particular, since the early nineties of the last 
century and in the last ten-fifteen years, a quiet revolution has taken place in the quantum 
arena. This revolution has progressively indicated and clarified that aspects once thought 
to be problematic, such as quantum nonseparability and "spooky" actions at a distance, are 
actually not only problems but rather some of the key ingredients that are allowing a deeper 
understanding of quantum mechanics, its applications to new and exciting fields of research 
(such as quantum information and quantum computation), and tremendous progress in the 
development of its mathematical and conceptual foundations. Among the key elements of the 
current re-foundation of quantum mechanics, entanglement certainly plays a very important 
role because it is a concept that can be mathematically qualified and quantified in a way that 
allows it to provide new and general characterizations of quantum properties, operations, and 
states. 

In the context of this special issue on Quantum Information, we will review the 
main aspects of entanglement in continuous variable systems, as they are currently 
understood. With the aim to be as self-contained as possible, we start with a tutorial 
summary of the structural properties of continuous variable systems -systems associated to 
infinite-dimensional Hilbert spaces whose possible applications in quantum information and 
communication tasks are gaining increasing interest and attention- focusing on the specially 
relevant family of Gaussian states. The reason why this review is mainly focused on the 
discussion of separability criteria and entanglement properties of Gaussian states is due to 
their great practical relevance in applications to quantum optics and quantum information, to 
the very clean framework that they allow for the study of the structure of nonlocal correlations, 
and to the obvious consequence that in the last years most studies and results on continuous- 
variable entanglement have been obtained for Gaussian states. 

In Sections 2 and 3 we give a self-contained introduction to phase-space and symplectic 
methods in the study of Gaussian states of infinite-dimensional bosonic systems, we discuss 
the covariance matrix formalism, and we provide a classification of pure and mixed Gaussian 
states according to the various forms that the associated covariance matrices can take. In 
section 4 we introduce and describe the separability problem, the inseparability criteria, and 
the machinery of bipartite entanglement evaluation and distillation in Gaussian states. In 
Section 5 we review some important results specific on two-mode Gaussian states, including 
the existence of extremally entangled states, minimal and maximal, at a given degree of 
mixedness, and the different orderings induced on the set of two-mode Gaussian states by 
different measures of entanglement such as the Gaussian entanglement of formation and 
the logarithmic negativity. In Section 6 we describe the unitary (and therefore reversible) 
localization of bipartite multimode entanglement to bipartite two-mode entanglement in 
fully symmetric and bisymmetric multimode Gaussian states, and the scaling of bipartite 
entanglement with the number of modes in general multimode Gaussian states. In Section 
7 we then discuss recent crucial advances in the understanding of entanglement sharing in 
multimode Gaussian states, including the proof of the monogamy inequality of distributed 
entanglement for all Gaussian states. Multipartite entanglement of Gaussian states is reviewed 
in Section 8 by discussing its qualification according to different classes of separability, and 
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the main consequences of the monogamy inequality, such as the quantification of genuine 
tripartite entanglement in three-mode Gaussian states via the residual Gaussian tangle, the 
promiscuous nature of entanglement sharing in Gaussian states with symmetry constraints, 
and the possible coexistence of unlimited bipartite and multipartite entanglement. The last two 
properties (promiscuity and coexistence of unbound bipartite and multipartite entanglement) 
are elucidated at length by discussing three- and four-mode Gaussian states endowed with 
full and partial symmetry constraints under the exchange of modes. Section 9 concludes this 
review with a brief discussion about recent advances in the qualification and quantification 
of entanglement in non Gaussian states, a field of research that is to a large extent yet to 
be fully explored, a summary on various applications of Gaussian entanglement in quantum 
information and computation, and on overview on open problems and current research 
directions. 



2. Introduction to continuous variable systems. 

A continuous variable (CV) system 0] O of N canonical bosonic modes is described 
by a Hilbert space H = (££) fc=1 Ti-k resulting from the tensor product structure of infinite 
dimensional Fock spaces 7ife's, each of them associated to a single mode. For instance, one 
can think of the non interacting quantized electromagnetic field, whose Hamiltonian describes 
a system of an arbitrary number N of harmonic oscillators of different frequencies, the modes 
of the field, 

N , 



H = ^hw k I a\a k + - J . (1) 

k=l ^ ' 

Here a k and at are the annihilation and creation operators of a photon in mode k (with 
frequency uiu), which satisfy the bosonic commutation relation 



a k ,a\, 



4>4' 



0. (2) 



Skk> , [fifcOfc'] 

From now on we will assume for convenience natural units with h = 2. The corresponding 
quadrature phase operators (position and momentum) for each mode are defined as 

qk = (Sfc + a\) , (3) 

Pk = {a k - a\)/i (4) 

We can group together the canonical operators in the vector 

R = (qi,Pi, ■ ■ ■ ,Qn,Pn) T , (5) 

which enables us to write in compact form the bosonic commutation relations between the 
quadrature phase operators, 

[R k ,Ri] =2in H , (6) 

where fl is the symplectic form 

« = -=(_°! I)- W 

fc=l ^ ' 

The space H y. is spanned by the Fock basis { | n) & } of eigenstates of the number operator 
fik = a\a>k> representing the Hamiltonian of the noninteracting mode via Eq. (|T|). The 
Hamiltonian of each mode is bounded from below, thus ensuring the stability of the system. 
For each mode k there exists a different vacuum state |0)/. S TLk such that afc|0)fe = 0. The 
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vacuum state of the global Hilbert space will be denoted by |0) = |Q)fe. In the single- 
mode Hilbert space Hk, the eigenstates of constitute the important set of coherent states 
H, which is overcomplete in Hk- Coherent states result from applying the single-mode Weyl 
displacement operator Dk to the vacuum |0)fc, \a)k = Dk(a)\0)k, where 

D k (a) = e a ^~ a * a * , (8) 

and the coherent amplitude a 6 (D satisfies &k\a)k = a\a)k- In terms of the Fock basis of 
mode k a coherent state reads 

00 n 

l«>*= ^ |Q|2 ]T^=|n) fc . < 9 ) 

n— 1 v 

Tensor products of coherent states for iV different modes are obtained by applying the N- 
mode Weyl operators to the global vacuum |0). For future convenience, we define the 
operators in terms of the canonical operators R, 

£>z = e^ Tf2 « , with £ G R 2N . (10) 
One then has |f) = D^O). 



2.1. Quantum phase-space picture. 

The states of a CV system are the set of positive trace-class operators {g} on the Hilbert space 
H = T~tk- However, the complete description of any quantum state g of such an infinite 

dimensional system can be provided by one of its s-ordered characteristic functions 

X s(0= Tr[^]e s "«H 2 / 2 , (11) 

with £ £ ]R 2Af , || . || standing for the Euclidean norm of H 2N . The vector £ belongs to the real 
27V-dimensional space T = (TR 2N , O), which is called phase space, in analogy with classical 
Hamiltonian dynamics. One can see from the definition of the characteristic functions that 
in the phase space picture, the tensor product structure is replaced by a direct sum structure, 
so that the iV-mode phase space T = r&, where Tk — (R 2 , u>) is the local phase space 
associated with mode k. 

The family of characteristic functions is in turn related, via complex Fourier transform, to 
the quasi-probability distributions W s , which constitute another set of complete descriptions 
of the quantum states 

W.{£) = \ ( k Xs {k) e lKTf2 « d 2N . (12) 

There exist states for which the function W s is not a regular probability distribution for 
any s, because it can be singular or assume negative values. Note that the value s = — 1 
corresponds to the Husimi 'Q-function' J6) W-i(£) — (C|p|0/ 7r an< ^ thus always yields 
a regular probability distribution. The case s = corresponds to the so called Wigner 
function Q, which will be denoted simply by W. Likewise, for the sake of simplicity, x 
will stand for the symmetrically ordered characteristic function xo- Finally, the case s = 1 
yields the singular P-representation, which was introduced, independently, by Glauber 1 8 1 and 
Sudarshan J9). 

The quasiprobability distributions of integer order W-%, Wo and Wi are respectively 
associated the antinormally ordered, symmetrically ordered and normally ordered expressions 
of operators. More precisely, if the operator O can be expressed as O = f(ak,a/ k ) for 
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Q 





Table 1. Schematic comparison between Hilbert-space and phase-space 
pictures for iV-mode continuous variable systems. 



k = 1,...,N, where / is a, say, symmetrically ordered function of the field operators, then 
onehaslfTOlfTn 

TT[g6] = f Wo(K)f( K )d 2N K, 

J-R 2N 

where f(n) = /(«;& + iKk+i,Kk — ink+i) and / takes the same form as the operatorial 
function previously introduced. The same relationship holds between W—\ and the 
antinormally ordered expressions of the operators and between W\ and the normal ordering. 
We also recall that the normally ordered function of a given operator is provided by its Wigner 
representation. This entails the following equalities for the trace 

l=Trf>= / W( K )d 2N K = X (0), (13) 
and for the purity lfl2l 

M = Tr, 2 =/ W\K)d™K=( |x(C)| 2 d 2 ^, (14) 

of a state g. These expressions will be useful in the following. 

The (symmetric) Wigner function can be written as follows in terms of the 
(unnormalized) eigenvectors \x) of the quadrature operators {qj} (for which q~j\x) — qj\x), 
x g TR, N . for 7 = 1 N) Ifl3l 

W{x,p) = \ [ (x- x'\g\x + x')e tx '- p d N x l , x,p€B. N (15) 

From an operational point of view, the Wigner function admits a clear interpretation in terms 
of homodyne measurements [14|: the marginal integral of the Wigner function over the 
variables p u . . . ,p N ,x x , . . . ,xjv-i 

/ W(x,p) d N pdxi . . . dx^v-i 

gives the probability of the results of homodyne detections on the remaining quadrature x n 

ma. 

Table [TJ summarizes the mapping of properties and tools between Hilbert and phase 
spaces. In the next section more properties and tools of phase space analysis will be introduced 
that are specially useful in the study of Gaussian states. 

3. Mathematical description of Gaussian states. 

The set of Gaussian states is, by definition, the set of states with Gaussian characteristic 
functions and quasi-probability distributions on the multimode quantum phase space. 
Gaussian states include, among others, coherent, squeezed, and thermal states. Therefore 
they are of central importance in quantum optics and in quantum information and quantum 
communication with CV systems [ 16 1. Their entanglement properties will thus be one of the 
main subjects of this review. 
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3.1. Covariance matrix formalism. 

From the definition it follows that a Gaussian state q is completely characterized by the 
first and second statistical moments of the quadrature field operators, which will be denoted, 

respectively, by the vector of first moments R = (j^Rx), (Ri), ■ ■ ■ , (Rn), (Rn)^j an d by the 
covariance matrix (CM) cr of elements 

a tj = ~ (RiRj + RjRi) - {R t ) (Rj) , . (16) 

First moments can be arbitrarily adjusted by local unitary operations, namely 
displacements in phase space, i.e. applications of the single-mode Weyl operator Eq. ([8]) to 
locally re-center the reduced Gaussian corresponding to each single mode (recall that the 
reduced state obtained from a Gaussian state by partial tracing over a subset of modes is still 
Gaussian). Such operations leave any informationally relevant property, such as entropy and 
entanglement, invariant. Therefore, from now on (unless otherwise stated) we will adjust first 
moments to without any loss of generality for the scopes of our analysis. 

With this position, the Wigner function of a Gaussian state can be written as follows in 
terms of phase-space quadrature variables 

wm —jBsr- <17) 

where R stands for the real phase-space vector (<Zi,Pi, ■ • ■ , q n ,Pn) G T. Therefore, despite 
the infinite dimension of the associated Hilbert space, the complete description of an arbitrary 
Gaussian state (up to local unitary operations) is given by the 2 N x 2 N CM cr. In the following 
cr will be assumed indifferently to denote the matrix of second moments of a Gaussian state, or 
the Gaussian state itself. In the language of statistical mechanics, the elements of the CM are 
the two-point truncated correlation functions between the 2 N canonical continuous variables. 
We notice also that the entries of the CM can be expressed as energies by multiplying them by 
the level spacing hujk, where to k is the frequency of each mode k. In this way Tr cr is related 
to the mean energy of the state, i.e. the average of the noninteracting Hamiltonian Eq. ((T). 

As the real cr contains the complete locally-invariant information on a Gaussian state, we 
can expect some constraints to exist to be obeyed by any bona fide CM, reflecting in particular 
the requirements of positive-semidefiniteness of the associated density matrix q. Indeed, such 
condition together with the canonical commutation relations imply 

a- + m>o, (18) 

Ineq. JT8l is the necessary and sufficient constraint the matrix cr has to fulfill to be a CM 
corresponding to a physical Gaussian state flTl [T8l . More in general, the previous condition 
is necessary for the CM of any, generally non Gaussian, CV state (characterized in principle 
by the moments of any order). We note that such a constraint implies cr > 0. Ineq. (fT8l is the 
expression of the uncertainty principle on the canonical operators in its strong, Robertson- 
Schrodinger form |[T9l |20l |2~T1 . 

For future convenience, let us define and write down the CM <t\...n of an iV-mode 
Gaussian state in terms of two by two submatrices as 

/ en £i,2 ••• £i,n \ 



cr 1...JV — 



4,2 



£n-i,n 

-l.N ' ' ' £ N-1,N 



(19) 
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Each diagonal block <r fc is respectively the local CM corresponding to the reduced state of 
mode k, for all k = 1, . . . , N, On the other hand, the off-diagonal matrices £j j encode the 
intermodal correlations (quantum and classical) between subsystems i and j. The matrices 
ei.j all vanish for a product state. 

In this preliminary overview, let us just mention an important instance of two-mode 
Gaussian state, the two-mode squeezed state \ip sq )i,j 



Ui,j(r) (|0)j($ |0)j) with squeezing 



factor rel, where the (phase-free) two-mode squeezing operator is given by 



Ui,j(r) = exp 



7 Vt~t - - \ 



(20) 



In the limit of infinite squeezing (r — ► oo), the state approaches the ideal Einstein-Podolsky- 
Rosen (EPR) state [22|, simultaneous eigenstate of total momentum and relative position of 
the two subsystems, which thus share infinite entanglement. The EPR state is unnormalizable 
and unphysical. However, in principle, an EPR state can be approximated with an arbitrarily 
high degree of accuracy by two-mode squeezed states with sufficiently large squeezing. 
Therefore, two-mode squeezed states are of key importance as entangled resources for 
practical implementations of CV quantum information protocols [1]. They thus play a central 
role in the subsequent study of the entanglement properties of Gaussian states. A two-mode 
squeezed state with squeezing r (also known in quantum optics as a twin-beam state) is 
described by a CM 

I cosh(2r) sinh(2r) \ 

cosh(2r) - sinh(2r) 

sinh(2r) cosh(2r) 

\ -sinh(2r) cosh(2r) / 



(r) 



(21) 



3.2. Symplectic operations. 

An important role in the theoretical and experimental manipulation of Gaussian states is 
played by unitary operations which preserve the Gaussian character of the states on which 
they act. These unitary Gaussian operations are all those generated by Hamiltonian terms 
at most quadratic in the field operators. As a consequence of the Stone- Von Neumann 
theorem, the so-called metaplectic representation entails that any such unitary operation at 
the Hilbert space level corresponds, in phase space, to a symplectic transformation, i.e. to a 
linear transformation 5 which preserves the symplectic form f2: 

s T ns = n . (22) 

Symplectic transformations on a 2iV-dimensional phase space form the (real) symplectic 
group Sp(2N,n)- Such transformations act linearly on first moments and by congruence on 
covariance matrices, er i— » ScrS T . Eq. (l22l implies Det5 = 1, VS G Sp( 2 N,B,)- Ideal 
beam splitters, phase shifters and squeezers are all described by some kind of symplectic 
transformation (see e.g. [15]). For instance, the two-mode squeezing operator Eq. (f20b 
corresponds to the symplectic transformation 



Sij(r) 



( 



coshr 


sinhr 






coshr 


— sinh r 



sinh r 


coshr 






— sinhr 


coshr 



(23) 



where the matrix is understood to act on the pair of modes i and j. In this way, the two-mode 
squeezed state, Eq. (|2TT >. can be obtained as o - *' (r) = Sij (r) ISjj (r) exploiting the fact that 
the CM of the two-mode vacuum state is the 4x4 identity matrix. 
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Another common unitary operation is the ideal (phase-free) beam splitter, whose action 
Bi j on a pair of modes i and j is defined as 



a,j — > a,; sm W — a., 



with a; being the annihilation operator of mode fc. A beam splitter with transmittivity r 
corresponds to a rotation of 9 = arccos y/r in phase space (6 = tt/4 corresponds to a balanced 
50:50 beam splitter, r = 1 /2), described by a symplectic transformation 



/ Vt o 7 o \ 

y/T v 7 !" 7 r 



(25) 



VT^T -V^ 

V V^~T -V? J 
Single-mode symplectic operations are easily introduced as linear combinations of planar 
(orthogonal) rotations and of single-mode squeezings of the form 

Sj(r)= diag(e r , e" r ) , (26) 

acting on mode j, for r > 0. 

In general, symplectic transformations in phase space are generated by the 
exponentiation of matrices written as JO, where J is antisymmetric [23 1. Such generators can 
be symmetric or antisymmetric. The operations Bij (t), Eq. d25l l, generated by antisymmetric 
operators are orthogonal and, acting by congruence on the CM cr, preserve the value of Tr cr. 
Since Tr cr gives the contribution of the second moments to the average of the Hamiltonian 
©fc ^k^k, these transformations are said to be passive (they belong to the compact subgroup 
of S'p(2Ar,R)). Instead, operations Sij(r), Eq. (l23l , generated by symmetric operators, are not 
orthogonal and do not preserve Tr cr (they belong to the non compact subgroup of 5p(2jv,R))- 
This mathematical difference between squeezers and phase-space rotations accounts, in a 
quite elegant way, for the difference between active (energy non preserving) and passive 
(energy preserving) optical transformations 1 241 . 

Let us remark that local symplectic operations belong to the group Sp(2, R)® . They 
correspond, at the Hilbert space level, to tensor products of unitary transformations, each 
acting on the state space of a single mode. It is useful to notice that the determinants of each 
2x2 submatrix of a A^-mode CM, Eq. (T% , are all invariants under local symplectic operations 
S G Sp(2, M) eAr . This mathematical property reflects the physical requirement that marginal 
informational properties and correlations between individual subsystems cannot be altered by 
local operations alone. 

3.2.1. Symplectic eigenvalues and invariants. A symplectic transformation of major 
importance is the one that diagonalizes a Gaussian state in the basis of normal modes. Through 
this decomposition, thanks to Williamson theorem ll25l . the CM of a A-mode Gaussian state 
can always be written in the so-called Williamson normal, or diagonal form 

cr = S T vS , (27) 
where S € Sp^N.R) an d v i s tne CM 

N 



v k 
v k 



(28) 



k=l 

corresponding to a tensor product state with a diagonal density matrix q® given by 
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where \n)f. denotes the number state of order n in the Fock space Hk- In the Williamson 
form, each mode with frequency ojk is a Gaussian state in thermal equilibrium at a temperature 
Tfe, characterized by a Bose-Einstein statistical distribution of the thermal photons rife, with 
average 

Uk — 1 1 
n k = = -r— -r . (30) 

2 expfe)-! 



s k B T k/ 

The N quantities u k 's form the symplectic spectrum of the CM er, and are invariant 
under the action of global symplectic transformations on the matrix er. The symplectic 
eigenvalues can be computed as the orthogonal eigenvalues of the matrix |if2<r| ll26l and 
are thus determined by N invariants of the characteristic polynomial of such a matrix ETl . 
One global symplectic invariant is simply the determinant of the CM (whose invariance is 
a consequence of the fact that Det S = 1 V5 G S'p(2JV,R))) which, once computed in the 
Williamson diagonal form, reads 

N 

Det<r = l[4. (31) 

fc=i 

Another important invariant under global symplectic operations is the so-called seralian A 
ll27ll . defined as the sum of the determinants of all 2 x 2 submatrices of a CM er, Eq. ( fT9l >. 
which can be readily computed in terms of its symplectic eigenvalues as 

JV 

A( CT )=^. (32) 

fc=l 

The invariance of A CT in the multimode case OTTl follows from its invariance in the case of 
two-mode states, proved in Ref. 11281 . and from the fact that any symplectic transformation 
can be decomposed as the product of two-mode transformations l29]| . 



3.2.2. Symplectic representation of the uncertainty principle. The symplectic eigenvalues 
Vk provide a powerful tool to access essential information on the properties of the Gaussian 
state er J28). For instance, let us consider the uncertainty relation dl~8b . Since the inverse 
of a symplectic operation is itself symplectic, one has from Eq. ( l22l . S~ lT ClS~ 1 = O, so 
that Ineq. ( fT8b is equivalent to v + ifl > 0. In terms of the symplectic eigenvalues v k the 
uncertainty relation then simply reads 

v k > 1 . (33) 

Inequality d33l is completely equivalent to the uncertainty relation ( fT8l provided that the CM 
er satisfies er > 0. 

Without loss of generality, one can rearrange the modes of a iV-mode state such that the 
corresponding symplectic eigenvalues are sorted in ascending order 

V— = V\ < V2 < ■ ■ • < Vn-i < v n = v+ . 

With this notation, the uncertainty relation reduces to v\ > 1. We remark that the full 
saturation of the uncertainty principle can only be achieved by pure TV-mode Gaussian states, 
for which 

Vi = 1 V*= l,...,N, 

meaning that the Williamson normal form of any pure Gaussian state is the vacuum |0) of 
the iV-mode Hilbert space H. Instead, mixed states such that fi<k = 1 and v%>k > 1> with 
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Table 2. Schematic comparison between Hilbert-space and phase-space pictures for A^-mode 
Gaussian states. The first two rows are taken from Table[T]and apply to general states of CV 
systems. The following rows are special to Gaussian states, relying on the covariance matrix 
description and the properties of the symplectic group. 



1 < k < N, only partially saturate the uncertainty principle, with partial saturation becoming 
weaker with decreasing k. Such states are minimum uncertainty mixed Gaussian states in 
the sense that the phase quadrature operators of the first k modes satisfy the Robertson- 
Schrodinger minimum uncertainty, while, for the remaining N — k modes, the state indeed 
contains some additional thermal correlations which are responsible for the global mixedness 
of the state. 

The symplectic rank H of a CM cr is defined as the number of its symplectic eigenvalues 
different from 1, corresponding to the number of non-vacua normal modes ||30l . A Gaussian 
state is pure if and only if H = 0. For mixed iV-mode states one has 1 < H < N 
depending on their degree of partial minimum-uncertainty saturation. This is in analogy 
with the standard rank of finite-dimensional (density) matrices, defined as the number of 
nonvanishing eigenvalues; in that case, only pure states q — nave rar, k 1, and mixed 

states have in general higher rank. As we will now show, all the informational properties of 
Gaussian states can be recast in terms of their symplectic spectra. 

A mnemonic summary of the main structural features of Gaussian states in the phase- 
space/CM description (including the definition of purity given in the next subsection) is 
provided by Table [2] 



3.3. Degree of information contained in a Gaussian state. 

3.3.1. Measures of information. The degree of information contained in a quantum state 
corresponds to the amount of knowledge that we possess a priori on predicting the outcome 
of any test performed on the state 1 3 1 1 . 

The simplest measure of such information is the purity of a quantum state q: 

Me)=Trg 2 . (34) 

For states belonging to a given Hilbert space H with dim7Y = D, the purity varies in the 
range 

1 

— < u < 1 . 

D ~ 

The minimum is reached by the totally random mixture; the upper bound is saturated by pure 
states. In the limit of CV systems (D — > oo), the minimum purity tends asymptotically 



Entanglement in continuous variable systems 



12 



to zero. Accordingly, the "impurity" or degree of mixedness of a quantum state g, which 
characterizes our ignorance before performing any quantum test on q, can be quantified by 
the functional 

S l (q) = ^ t (1-^ = ^ t {1-Tt Q 2 ) . (35) 

The linear entropy Sl (ranging between and 1) defined by Eq. d35l) is a very useful 
measure of mixedness in quantum mechanics and quantum information theory due to its direct 
connection with the purity and its computational simplicity. 

In general, the degree of mixedness of a quantum state q can be characterized completely 
by the knowledge of all the associated Schatten p-norms 11321 

\\ Q \\ p = (Tr\g\ p )* = (Tre p )p , withp > 1. (36) 

In particular, the case p = 2 is directly related to the purity p,, Eq. d34| >, as it is 
essentially equivalent (up to a normalization) to the linear entropy, Eq. d35l >. The p-norms 
are multiplicative on tensor product states and thus determine a family of non-extensive 
"generalized entropies" S p ll33l[34l . defined as 

S p ( Q )= 1 ~ Tr / , p>l. (37) 

p- 1 

The generalized entropies S p 's range from for pure states to l/(p— 1) for completely mixed 
states with fully degenerate eigenspectra. We also note that, for any given quantum state, 
S p is a monotonically decreasing function of p. Finally, another important class of entropic 
measures includes the Renyi entropies [|35l 

S» = ^^, P>1- (38) 
F 1 — p 

It can be shown that 



limS p (g) = lim S*(g) = -Tr(glogg) . (39) 

Therefore, even the von Neumann entropy 

S v (g) = -Tv g log g (40) 

can be defined, in terms of p-norms, as a suitable limit within the framework of generalized 
entropies. 

The von Neumann entropy is subadditive ll36l . Consider a bipartite system S (described 
by the Hilbert space Ti = Hi (8 H-2) in the state g. Then 

Sv(g) <S v (gi) + S v (g 2 ), (41) 

where gi t 2 are the reduced density matrices g\ t 2 — Tr2.i g associated to subsystems S\^- For 
states of the form g® = gi® g 2 , Eq. d4T1 ) is saturated, yielding that the von Neumann entropy 
is additive on tensor product states: 

Sv(ei<8ea) = S v ( gi ) + S v (g 2 ) . (42) 

The purity, Eq. d34l l. is instead multiplicative on product states, as the trace of a product 
equates the product of the traces: 

h(qi ® ei) = Kei) ■ m(^) ■ (43) 

Considering a joint classical probability distribution over the variables X and Y, one has 
for the Shannon entropy, 

S(X,Y) >S(X), S(Y). (44) 
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The Shannon entropy of a joint probability distribution is always greater than the Shannon 
entropy of each marginal probability distribution, meaning that there is less information in 
a global classical system than in any of its parts. On the other hand, consider a bipartite 
quantum system in a pure state q — ■ We have then for the von Neumann entropies: 

Sv(q) = 0, Sv(qi) = Sv{Q2) > 0. The global state g has been prepared in a well defined 
way, but if we measure local observables on the subsystems, the measurement outcomes are 
unavoidably affected by random noise. One cannot reconstruct the whole information about 
how the global system was prepared in the state g (apart from the trivial instance of g being 
a product state g = gi £2), by only looking separately at the two subsystems. Information 
is rather encoded in nonlocal and non factorizable quantum correlations - entanglement - 
between the two subsystems. This clearly highlights the fundamental difference between 
classical and quantum distributions of information. 



3.3.2. Entropic measures for Gaussian states. We now illustrate how to evaluate the 
measures of information (or lack thereof) defined above for Gaussian states (37). 

The generalized purities Tr g p defined by Eq. d36l l are invariant under global unitary 
operations. Therefore, for any TV-mode Gaussian state they are only functions of the 
symplectic eigenvalues Vf. of cr. In fact, a symplectic transformation acting on a is embodied 
by a unitary (trace preserving) operator acting on g, so that Tr g p can be easily computed on 
the Williamson diagonal state v of Eq. d28| >. One obtains 11371 

N 

K Q" = l[9p(vi) , (45) 

1=1 

where 

( 2P 
9p[X) ~ (x + 1)p -{x- 1)p • 

A remarkable consequence of Eq. ( l43T > is that 

^ = T^rvm^- (46) 

Regardless of the number of modes, the purity of a Gaussian state is fully determined by the 
global symplectic invariant Det cr alone, Eq. PTT l. We recall that the purity is related to the 
linear entropy Sl via Eq. (1351 . which in CV systems simply becomes <Sx = 1 — A second 
consequence of Eq. ( l45l l is that, together with Eqs. (T37T > and (f39b . it allows for the computation 
of the von Neumann entropy Sy, Eq. (|40| |. of a Gaussian state g, yielding [28 1 



where 



N 

S v ( g ) = J2f(vi) , (47) 

i=l 

fix) = — log ^— J -—log [— ) ■ (48) 

Such an expression for the von Neumann entropy of a Gaussian state was first explicitly given 
in Ref. 11381 . Notice that Sy diverges on infinitely mixed CV states, while Sl is normalized 
to 1. Let us remark that, clearly, the symplectic spectrum of single-mode Gaussian states, 
which consists of only one eigenvalue v\, is fully determined by the invariant Det cr = v\. 
Therefore, all the entropies S p 's (and Sv as well) are just increasing functions of Det cr (i.e. of 
Sl) and induce the same hierarchy of mixedness on the set of one-mode Gaussian states. This 
is no longer true for multi-mode states, even in the simplest instance of two-mode states l37l . 
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3.4. Standard forms of special Gaussian states. 

The symplectic analysis applied on Gaussian CMs implies that some of them can be suitably 
reduced under local operations. Such reductions go under the name of standard forms. 



3.4.1. Pure states: Phase-space Schmidt decomposition. In general, pure Gaussian states of 
a bipartite CV system admit a physically insightful decomposition at the CM level 139, 401 I4TI . 
which can be regarded as the direct analogue of the Schmidt decomposition for pure discrete- 
variable states [42|. Let us recall what happens in the finite-dimensional case. With respect 
to a bipartition of a pure state \iP)a\b into two subsystems Sa and Sb, one can diagonalize 
(via an operation Ua <8> Ub which is local unitary according to the considered bipartition) 
the two reduced density matrices qa,b, to find that they have the same rank and exactly the 
same nonzero eigenvalues {\k} (k = 1 . . . min{dim7i^4, dim7is}). The reduced state of 
the higher-dimensional subsystem (say Sb) will accomodate (dim?4 — dim 7^) additional 

O's in its spectrum. The state \iP)a\b takes thus the Schmidt form \ip) = Ylk=i ^fcl u fc> v k)- 

Looking at the mapping provided by Table [2]one can deduce what happens for Gaussian 
states. Given a Gaussian CM rr a\b of an arbitrary number N of modes, where subsystem Sa 
comprises Na modes and subsystem Sb Nb modes (with Na + Nb = N), then one can 
perform the Williamson decomposition Eq. ( f27T > in both reduced CMs (via a local symplectic 
operation Sa © Sb), to find that they have the same symplectic rank, and the same non-unit 
symplectic eigenvalues {vk} (k = 1 . . . minjA^, Nb})- The reduced state of the higher- 
dimensional subsystem (say Sb) will accomodate (Nb — Na) additional l's in its symplectic 
spectrum. With respect to an arbitrary A\B bipartition, therefore, the CM cr p of any pure 
TV-mode Gaussian state is locally equivalent to the form <r p s — (Sa © Sb)c p (Sa © Sb) 1 , 
with 



r P - 



V 



-«4 



(49) 



Here each element denotes a 2 x 2 submatrix, in particular the diamonds (o) correspond to 
null matrices, 1 to the identity matrix, and 



















1 



The matrices Ck contain the symplectic eigenvalues i>k ^ 1 of both reduced CMs. By 
expressing them in terms of hyperbolic functions, Vk = cosh(2rfc) and by comparison with 
Eq. (f2Tb . one finds that each two-mode CM 

C fe S k 
S k Ck 

encoding correlations between a single mode from Sa and a single mode from Sb, is a two- 
mode squeezed state with squeezing r^. Therefore, the Schmidt form of a pure iV-mode 
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Gaussian state with respect to a (Na x iV^-mode bipartition (with Nq > Na) is that of a 
direct sum II391I4T1 

N A N 

^ = 0<^(n) <r° k , (50) 

i=l fc=2AT A +l 

where mode i S 5a, mode j = i + Na € Sb, and er° = I2 is the CM of the vacuum 
state of mode k 6 5s . This corresponds, on the Hilbert space level, to the product 
of two-mode squeezed states, tensor additional uncorrelated vacuum modes in the higher- 
dimensional subsystem (Sb in our notation) Boll . The phase-space Schmidt decomposition 
is a very useful tool both for the understanding of the structural features of Gaussian states 
in the CM formalism, and for the evaluation of their entanglement properties. Notice that 
the validity of such a decomposition can be extended to mixed states with a fully degenerate 
symplectic spectrum, i.e. a Williamson normal form proportional to the identity BOl I4T1 . As 
a straightforward consequence of Eq. d50l ). any pure two-mode Gaussian state is equivalent, 
up to local unitary operations, to a two-mode squeezed state of the form Eq. (f2TT> . 

We will now show that for (generally mixed) Gaussian states with some local symmetry 
constraints, a similar phase-space reduction is available, such that multimode properties (like 
entanglement) can be unitarily reduced to two-mode ones. 



3.4.2. Symmetric and bisymmetric states. Very often in quantum information, and in 
particular in the theory of entanglement, peculiar roles are played by symmetric states, that 
is, states that are either invariant under a particular group of transformations — like Werner 
states of quefits [43 1 — or under permutation of two or more parties in a multipartite system, 
like ground and thermal states of translationally invariant Hamiltonians (e.g. of harmonic 
lattices) l44]| . Here we will introduce classes of Gaussian states globally invariant under all 
permutation of the modes (fully symmetric states), or locally invariant in each ofvthe two 
subsystems across a global bipartition of the modes (bisymmetric states). The properties of 
their symplectic spectrum allow to determine explicitely a standard form in both cases. Here 
we will limit ourself to list and collect the results that will be useful for the computation and 
exploitation of entanglement in the corresponding states. The detailed discussions and all the 
rigorous proofs can be found in Ref. B31 . Unless explicitly stated, in the following we will 
be dealing with generally mixed states. 

A multimode Gaussian state q is "fully symmetric" if it is invariant under the exchange 
of any two modes. In the following, we will consider the fully symmetric M-mode and TV- 
mode Gaussian states g a M and Qqn, with CMs <j a M and ct^n . Due to symmetry, we have 
that 

/ a e ■■■ e \ ( (3 C ••• C \ 



a 



\ 



£ 

a 



T fi- 



P 

c 



(51) 



■ c 

V c ••• c p J 

where a, e, /3 and £ are 2 x 2 real symmetric submatrices (the symmetry of e and £ stems 
again from the symmetry under the exchange of any two modes). 

A fully symmetric iV-mode Gaussian states with CM ct^n admits the following standard 
form [46 , 45]. The 2x2 blocks (3 and £ of ct^n, defined by Eq. dBTl l. can be brought by means 
of local, single-mode symplectic operations S € Sp®^^ into the form (3 = diag (b, b) and 
C = diag (zi, Z2). Obviously, analogous results hold for the M-mode CM cr a M of Eq. ( IBTT l. 
whose 2x2 submatrices can be brought to the form a = diag (a, a) and e = diag (ei , e 2 ). 
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Let us now generalize this analysis to the (M + TV) -mode Gaussian states with CM er, 
which results from a correlated combination of the fully symmetric blocks cr a M and ct^n : 

T n M r 



T T er 



(52) 



where T is a 2M x 2N real matrix formed by identical 2x2 blocks 7. Clearly, T is 
responsible for the correlations existing between the A/-mode and the iV-mode parties. Once 
again, the identity of the submatrices 7 is a consequence of the local invariance under mode 
exchange, internal to the A/-mode and iV-mode parties. States of the form of Eq. ( |52l will be 
henceforth referred to as "bisymmetric" ll46l l45l . By evaluating the symplectic spectrum of 
such a bisymmetric state, and by comparing it with the spectra of the reduced fully symmetric 
CMs, one can show the following central result [45 1, which applies to all (generally mixed) 
bisymmetric Gaussian states, and is somehow analogous to — but independent of — the 
phase-space Schmidt decomposition of pure Gaussian states (and of mixed states with fully 
degenerate symplectic spectrum). The bisymmetric (M + N)-mode Gaussian state with CM 
cr, Eq. ( 1521 ) can be brought, by means of a local unitary (symplectic) operation with respect 
to the M X N bipartition with reduced CMs <r a M and ct^n, to a tensor product of single- 
mode uncorrelated states and of a two-mode Gaussian state comprised of one mode from the 
M-mode block and one mode from the N -mode block. 

We will explore the consequences of such a unitary localization on the entanglement 
properties of bisymmetric states, in the next sections. Let us just note that fully symmetric 
Gaussian states, Eq. ( BTT l. are special instances of bisymmetric states with respect to any global 
bipartition of the modes. 



4. Theory of bipartite entanglement for Gaussian states. 

4.1. Entanglement and nonlocality. 

When two quantum systems come to interact, some quantum correlation is established 
between the two of them. This correlation persists even when the interaction is switched 
off and the two systems are spatially separatee^ If we measure a local observable on the 
first system, its state collapses in an eigenstate of that observable. In the ideal case of no 
environmental decoherence, the state of the second system is instantly modified. The question 
then arises about what is the mechanism responsible for this "spooky action at a distance" 

(22. 

Suppose we have a bipartite or multipartite quantum state. An apparently innocent 
question like 

Does this state contain quantum correlations? 

turns out to be very hard to answer II481I491I501 , The first step towards a solution concerns a 
basic understanding of what such a question really means. 

One may argue that a system contains quantum correlations if the observables associated 
to the different subsystems are correlated, and their correlations cannot be reproduced with 
purely classical means. This implies that some form of inseparability or non factorizability is 
necessary to properly take into account those correlations. For what concerns globally pure 
states of the composite quantum system, it is relatively easy to control if the correlations are 
of genuine quantum nature. In particular, it is enough to check if a Bell-CHSH inequality 

§ Entanglement can be also created without direct interaction between the subsystems, via the common interaction 
with other assisting parties . This and related mechanisms define the so-called entanglement swapping | 47 1. 
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I15T1I521 is violated [5311 . to conclude that a pure quantum state is entangled. There are in fact 
many different criteria to characterize entanglement, but all of them are practically based on 
equivalent forms of nonlocality in pure quantum states. 

This scheme fails with mixed states. Contrary to a pure state, a statistical mixture 
can be prepared in (generally infinitely) many different ways. Being it impossible to 
reconstruct the original preparation of the state, one cannot extract all the information it 
contains. Accordingly, there is no completely general and practical criterion to decide whether 
correlations in a mixed quantum state are of classical or quantum nature. Moreover, different 
manifestations of quantum inseparability are in general not equivalent. For instance, one pays 
more (in units of singlets ( 1 00) + 1 1 1) )/ \/2) to create an entangled mixed state g (entanglement 
cost 1 54 1) than to reconvert g into a product of singlets (distillable entanglement |54|) via 
local operations and classical communication (LOCC) 1531 . An important example has been 
introduced by Werner 11431 . who defined a parametric family of mixed states (known as Werner 
states) which, in some range of the parameters, are entangled (inseparable) without violating 
any Bell inequality on local realistic models. These states thus admit a description in terms of 
local hidden variables. It is indeed an open question in quantum information theory to prove 
whether any entangled state violates some kind of Bell-type inequality 15611571 . 

In fact, entanglement and nonlocality are different resources from an operational point 
of view (58). This can be understood within the general framework of no-signalling theories 
which exhibit stronger nonlocal features than those of quantum mechanics. Let us briefly 
recall what is intended by nonlocality according to Bell (59): there exists in Nature a channel 
that allows one to distribute correlations between distant observers, such that the correlations 
are not already established at the source, and the correlated random variables can be created 
in a configuration of space-like separation, i.e. no normal signal (assuming no superluminal 
transmission) can be the cause of the correlations ll22l I5T1 . A convenient understanding of 
the phenomenon of nonlocality is then given by quantum mechanics, which describes the 
channel as a pair of entangled particles. But such interpretation is not the only possible one. 
In recent years, there has been a growing interest in constructing alternative descriptions of 
this channel, mainly assuming a form of communication [60 1, or the use of an hypothetical 
"nonlocal machine" [61 1 able to violate the CHSH inequality [52 1 up to its algebraic value 
of 4 (while the local realism threshold is 2 and the maximal violation admitted by quantum 
mechanics is 2y2, the Cirel'son bound [62 1). Looking at these alternative descriptions allows 
to quantify how powerful quantum mechanics is by comparing its performance with that of 
other resources [ 63 1 . 

4.2. Entanglement of pure states. 

It is now well understood that entanglement in a pure bipartite quantum state is equivalent 
to the degree of mixedness of each subsystem. Accordingly, it can be properly quantified by 
the entropy of entanglement Ev(\ip)), defined as the von Neumann entropy, Eq. d40l ). of the 
reduced density matrices l64l . 

d 

E v {\i>)) = Sv(gi) = S v (g 2 ) = -J2 A' logA 2 fe . (53) 

fc=i 

The entropy of entanglement is by definition invariant under local unitary operations 

^((fliStfaM) =E v (\i>)), (54) 

and it can be shown [65 1 that Ev(\ip)) cannot increase under LOCC performed on the state 
this is a fundamental physical requirement as it reflects the fact that entanglement cannot 
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be created via LOCC only [66, 67 1. It can be formalized as follows. Let us suppose, starting 
with a state \t/j) of the global system S, to perform local measurements on Si and 6> 2 > an d 
to obtain, after the measurement, the state \ipi) with probability pi, the state \(p 2 ) with 
probability and so on. Then 

Ev(\1>))>52PkEv(\<Pk)). (55) 

k 

Note that entanglement cannot increase on average, that is nothing prevents, for a given 
k, that Ev(\fk)) > Ev(\ip))- The concept of entanglement distillation is based on this 
fact |68, 54l l69ll : with a probability pu, it is possible to increase entanglement via LOCC 
manipulations. 

4.3. Entanglement of mixed states. 

A mixed state g can be decomposed as a convex combination of pure states, 

Q = ^2pk\ipk)(ipk\ ■ (56) 
k 

Eq. d56l ) tells us how to create the state described by the density matrix g: we have to prepare 
the state \ipi) with probability p\, the state with probability pa> etc. For instance, we 
could collect N copies (TV ^> 1) of the system, prepare n% ~ Npk of them in the state \ipk), 
and pick a copy at random. 

The difficulty lies in the fact that the decomposition of Eq. ( TSoT ) is not unique: apart from 
pure states, there exist infinitely many decompositions of a generic g, meaning that the mixed 
state can be prepared in infinitely many different ways. This has important consequences 
on the determination of mixed-state entanglement. Suppose that for a bipartite system in 
a mixed state we detect, by local measurements, the presence of correlations between the 
two subsystems. Given the ambiguity on the state preparation, we cannot know a priori 
if those correlations arose from a quantum interaction between the subsystems (meaning 
entanglement) or were induced by means of LOCC (meaning classical correlations). It is 
thus clear that a mixed state can be defined separable (classically correlated) if there exists 
at least one way of engineering it by LOCC; on the other hand it is entangled (quantumly 
correlated) if, among the infinite possible procedures for its preparation, there is no one which 
relies on LOCC alone l43l . 

However,deciding separability according to the above definition would imply checking 
all the infinitely many decompositions of a state g and looking for at least one, expressed 
as a convex combination of product states, to conclude that the state is not entangled. This 
is clearly impractical. For this reason, several operational criteria have been developed in 
order to detect entanglement in mixed quantum states 170117111721 . Some of them, of special 
relevance to Gaussian states of CV systems, are discussed in the following. 

4.4. Separability and distillability of Gaussian states. 

4.4.1. PPT criterion. One of the most powerful results to date in the context of separability 
criteria is the Peres-Horodecki condition f73l l74l . It is based on the operation of partial 
transposition of the density matrix of a bipartite system, obtained by performing transposition 
with respect to the degrees of freedom of one subsystem only. Peres criterion states that, if a 
state g s is separable, then its partial transpose gj 1 (with respect e.g. to subsystem Si) is a valid 
density matrix, in particular positive semidefinite, gj 1 > 0. Obviously, the same holds for gj 2 . 
Positivity of the partial transpose (PPT) is therefore a necessary condition for separability 
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||73ll . The converse (i.e. g Tl > g separable) is in general false, but is has been proven 
true for low-dimensional systems, specifically bipartite systems with Hilbert state space of 
dimensionality 2x2 and 2 x 3. In these cases the PPT property is equivalent to separability 
ll74ll . For higher dimensional tensor product structures of Hilbert spaces, PPT entangled states 
(with g Tl > 0) have been shown to exist ll75l . These states are known as bound entangled 
IT76ll , because their entanglement cannot be distilled to obtain maximally entangled states. The 
existence of bound entangled (undistillable) states with negative partial transposition has been 
conjectured [77, 78 1, but at present there is not yet evidence of this property ||56l . 

Recently, the PPT criterion has been generalized to continuous variable systems by 
Simon lfl3ll . who showed how the operation of transposition acquires in infinite-dimensional 
Hilbert spaces an elegant geometric interpretation in terms of time inversion (mirror reflection 
of the momentum operator). The PPT criterion is necessary and sufficient for the separability 
of all (1 x iV)-mode Gaussian states lfl3l l79l [80l . Under partial transposition, the CM (Ta\b, 
where subsystem Sa groups N A modes, and subsystem Sb is formed by Nb modes, is 
transformed into a new matrix 

& a\b = 9a\b & a\b 9a\b , (57) 

with 

9 A]B = diag{l, -1, 1, -1, ... , 1, -1, 1,1,1,1,..., 1, 1} . 

V v ' V v ' 

2N A 2N B 

Referring to the notation of Eq. (1191 . the partially transposed matrix <j a\b differs from er a\b 
by a sign flip in the determinants of the intermodal correlation matrices, Det £y ■, with modes 
i G Sa and modes j € sb- 

The PPT criterion yields that a Gaussian state a a\b (with Na — 1 and Nb arbitrary) is 
separable if and only if the partially transposed cr a\b is a bona fide CM, that is it satisfies the 
uncertainty principle Eq. ( fT8l . 

a A[B +itt>Q. (58) 

This property in turn reflects the positivity of the partially transposed density matrix g TA 
associated to the state g. For Gaussian states with Na > 1, not endowed with special 
symmetry constraints, the PPT condition is only necessary for separability, as bound entangled 
Gaussian states, whose entanglement is undistillable, exist already in the instance Na = 
Nr = 2 Il80l. 

We have demonstrated the existence of "bisymmetric" (Na + -/Vs)-mode Gaussian states 
for which PPT is again equivalent to separability [45 1 . In view of the invariance of the PPT 
criterion under local unitary transformations, and considering the results of section 13.4.21 on 
the unitary localization of bisymmetric Gaussian states, it is immediate to verity that the 
following property holds B31 . For generic Na X Ns-mode bipartitions, the positivity of 
the partial transpose (PPT) is a necessary and sufficient condition for the separability of 
bisymmetric (Na + Ng)-mode mixed Gaussian states of the form Eq. \52h In the case of 
fully symmetric mixed Gaussian states, Eq. ([51$ , of an arbitrary number of bosonic modes, 
PPT is equivalent to separability across any global bipartition of the modes. 

This statement generalizes to multimode bipartitions the equivalence between 
separability and PPT for 1 x N bipartite Gaussian states [ 80 1 . In particular, it implies that 
no bisymmetric bound entangled Gaussian states can exist ifSQl [8~T1 and all the Na x Nb 
multimode block entanglement of such states is distillable. Moreover, it justifies the use of 
the negativity and the logarithmic negativity as measures of entanglement for this class of 
multimode Gaussian states [6] 
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Physical 


Separable 


density matrix 


Q > 


Q Ta > 


covariance matrix 


cr + itt > 


fr + itt > 


symplectic spectrum 


Vh > 1 


ffe > 1 



Table 3. Schematic comparison between the conditions of existence and the conditions of 
separability for Gaussian states, as expressed in different representations. The second column 
qualifies the PPT condition, which is always implied by separability, and equivalent to it in 
1 X N, and in M X N bisymmetric Gaussian states. 



The PPT criterion has an elegant symplectic representation. The partially transposed 
matrix cr of any TV-mode Gaussian CM is still a positive and symmetric matrix. As such, it 
admits a Williamson normal-mode decomposition |25|, Eq. d27l ), of the form 

e> = S T v>S , (59) 

where S £ Sp^N.m.) and & is tne CM 

'-©(»*)■ 

k—l x ' 

The N quantities tVs are the symplectic eigenvalues of the partially transposed CM cr. 
The symplectic spectrum {vu} of cr encodes the structural and informational properties of a 
Gaussian state. The partially transposed spectrum {vk} encodes the qualitative (and, to some 
extent, quantitative — see section |43I ) characterization of entanglement in the state. The PPT 
condition d58l ), i.e. the uncertainty relation for cr, can be equivalently recast in terms of the 
parameters j>/.'s as 

h > 1 • (61) 

We can, without loss of generality, rearrange the modes of a TV-mode state such that the 
corresponding symplectic eigenvalues of the partial transpose cr are sorted in ascending order 

V- = V\ < 2>2 < ■ • • < VN-\ < Vfq = v + , 

in analogy to what done in section [3.2.21 for the spectrum of cr. With this notation, the PPT 
criterion across an arbitrary bipartition reduces to i>\ > 1 for all separable Gaussian states. If 
v\ < 1, the corresponding Gaussian state cr is entangled. The symplectic characterization of 
physical and PPT Gaussian states is summarized in Table|3] 

The distillability problem for Gaussian states has been solved in quite general terms [ 8 1 1 : 
the entanglement of any non-PPT bipartite Gaussian state is distillable by LOCC. However, 
we remind that this entanglement can be distilled only resorting to non Gaussian LOCC |82|, 
since distilling Gaussian states with Gaussian operations is impossible [83, 84l l85l . 



4.4.2. Additional criteria for separability. Let us briefly mention alternative characteriza- 
tions of separability for Gaussian states. 

For a general Gaussian state of any Na x Nb bipartition, a necessary and sufficient 
condition states that a CM cr corresponds to a separable state if and only if there exists a 
pair of CMs a a and r/ B , relative to the subsystems Sa and sb respectively, such that the 
following inequality holds ll80l . cr > ct a © ctb- This criterion is not very useful in practice. 
Alternatively, one can introduce an operational criterion based on iterative applications of a 
nonlinear map, that is independent of (and strictly stronger than) the PPT condition [86|, and 
completely qualifies separability for all bipartite Gaussian states. 
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Another powerful tool establish the separability of quantum states is given by the so- 
called entanglement witnesses. A state q is entangled if and only if there exists a Hermitian 
operator W such that Tr VV g < and Tr VV a > for any state a € V, where V C H is 
the convex and compact subset of separable states ll74l l57l . The operator W is the witness 
responsible for detecting entanglement in the state g. According to the Hahn-Banach theorem, 
given a convex and compact set T> and given g ^ V, there exists an hyperplane which 
separates g from T>. Optimal entanglement witnesses induce an hyperplane which is tangent 
to the set D (87). A sharper detection of separability can be achieved by means of nonlinear 
entanglement witnesses, curved towards the set T> l88l . A comprehensive characterization 
of linear and nonlinear entanglement witnesses is available for CV systems |89l , and can be 
efficiently applied to the detection of separability in arbitrary Gaussian states, both in the 
bipartite and in the multipartite context. 

Finally, several operational criteria have been developed, that are specially useful in 
experimental settings. They are based on the violation of inequalities involving combinations 
of variances of canonical operators, and their validity ranges from the two-mode [79 1 to the 
multimode setting [90 1 . 

4.5. Quantification of bipartite entanglement in Gaussian states. 

The question of the quantification of bipartite entanglement for general (pure and mixed) 
cannot be considered completely solved yet. We have witnessed a proliferation of 
entanglement measures, each one motivated by specific contexts in which quantum 
correlations play a central role, and accounting for different, and in some cases inequivalent, 
operational characterizations and orderings of entangled states. Detailed expositions on the 
subject can be found e.g. in Refs. fTO |9T] l50l l92l . 

The first natural generalization of the quantification of entanglement to mixed states is 
certainly the entanglement of formation Ep(g) J54), defined as the convex-roof extension |93| 
of the entropy of entanglement Eq. 653\ . i.e. the weighted average of pure-state entanglement, 



minimized over all decompositions of the mixed state g = J2kPk\i J k)(4'k\- This is clearly 
an optimization problem of formidable difficulty, and an explicit solution is known only for 
the mixed states of two qubits J94'|, and for highly symmetric states like Werner states and 
isotropic states in arbitrary dimension [95 , 96 1. In CV systems, an explicit expression for the 
entanglement of formation is available only for symmetric, two-mode Gaussian states [97|. 
To date, the additivity of the entanglement of formation remains an open problem [56|. 

4.5.7. Negativities. An important class of entanglement monotones is defined by the 
negativities, which quantify the violation of the PPT criterion for separability (see section 
14.4. U . i.e. how much the partial transposition of q fails to be positive. The negativity N(g) 
|98]|99) is defined as 




(62) 



M{g) 




(63) 



where 



TrVotO 



(64) 
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is the trace norm of the operator O. The negativity has the advantage of being a computable 
measure of entanglement, being 



N{q) =max|0,-^A-| , (65) 

where the {A^~ }'s are the negative eigenvalues of the partial transpose. 

The negativity can be defined for CV systems as well [100], even though a related 
measure is more often used, the logarithmic negativity Ej^{q) H100ll99l . 

Em(q) = log ||e T ' ||i = iog[i + 2AT( e )] . (66) 

The logarithmic negativity is additive and, despite not being convex, is a full entanglement 
monotone under LOCC 11011 : it is an upper bound for the distillable entanglement |99|, 
Ej^{q) > Ed(q), and coincides with the entanglement cost under operations preserving the 
positivity of the partial transpose [ 1 02 1 . Both the negativity and the logarithmic negativity fail 
to be continuous in trace norm on infinite-dimensional Hilbert spaces; however, this problem 
can be circumvented by restricting to physical states of finite mean energy [ 103 1. 

The great advantage of the negativities is that they are easily computable for general 
Gaussian states; they provide a proper quantification of entanglement in particular for 
arbitrary 1 X N and bisymmetric M X N Gaussian states, directly quantifying the degree 
of violation of the necessary and sufficient PPT criterion for separability, Eq. doTb . Following 
IU00l|2"6l[3"7ll46l , the negativity of a Gaussian state with CM cr is given by 

f idlfcP* 1 -!). iovk:v k <l. 
M{<t) = { (67) 
[ if v, > 1 Mi . 

Here the set {vk] is constituted by the symplectic eigenvalues of the partially transposed CM 
er. Accordingly, the logarithmic negativity reads 

{- Efc lo g^ for fc : i> fc < 1 . 
(68) 
if £>i > 1 Vi . 

The following lemma is quite useful for the interpretation and the computation of the 
negativities in Gaussian states [21 ]. In a [Na + -/Vs)-mode Gaussian state with CM a 'a\b< at 
most 

N min = rmn{N A , N B } (69) 

symplectic eigenvalues of the partial transpose cr a\b can violate the PPT inequality d6TT > 
with respect to a Na x Nb bipartition. Thanks to this result, in all 1 x N Gaussian states and 
in all bisymmetric M x N Gaussian states (whose symplectic spectra exhibit degeneracy, see 
section 13.4. 21 i. the negativities are quantified in terms of the smallest symplectic eigenvalue 
V- of the partially transposed CM alone. For > 1 the state is separable, otherwise it is 
entangled; the smaller the more entangled is the corresponding Gaussian state. In the 
limit of a vanishing partially transposed symplectic eigenvalue, i>- — > 0, the negativities 
grow unboundedly. In the special instance of two-mode Gaussian states, such a result had 
been originally derived in II 1 041 [371 . 
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4.5.2. Gaussian convex-roof extended measures. It is possible to define a family of 
entanglement measures exclusively defined for Gaussian states of CV systems [105|. The 
formalism of Gaussian entanglement measures (Gaussian EMs) has been introduced in 
Ref. [106 1 where the "Gaussian entanglement of formation" has been defined and analyzed. 
The framework developed in Ref. 01061 is general and allows to define generic Gaussian 
EMs of bipartite entanglement by applying the Gaussian convex roof, that is, the convex roof 
over pure Gaussian decompositions only, to any bona fide measure of bipartite entanglement 
defined for pure Gaussian states. 

The original motivation for the introduction of Gaussian EMs stems from the fact that 
the entanglement of formation J54), defined by Eq. (|62| >. involves a nontrivial minimization 
of the pure-state entropy of entanglement over convex decompositions of bipartite mixed 
Gaussian states in ensemble of pure states. These pure states may be, in principle, non 
Gaussian states of CV systems, thus rendering the analytical solution of the optimization 
problem in Eq. (l62b extremely difficult even in the simplest instance of one mode per party. 
Nevertheless, in the special subset of two-mode symmetric mixed Gaussian states, the optimal 
convex decomposition of Eq. d62l has been exactly determined, and is realized in terms of 
pure Gaussian states [97 1 . Apart from that case (which will be discussed in section l5~.2.21 i. the 
determination of the entanglement of formation of nonsymmetric two-mode Gaussian states 
(and more general Gaussian states) is an open problem in the theory of entanglement |56|. 
However, inspired by the results achieved on two-mode symmetric states, one can at first 
try to restrict the problem only to decompositions into pure Gaussian states. The resulting 
measure, named as Gaussian entanglement of formation in Ref. II 1 061 . is an upper bound to 
the true entanglement of formation and coincides with it for symmetric two-mode Gaussian 
states. 

In general, we can define a Gaussian EM Gb as follows. For any pure Gaussian state ip 
with CM er p , one has 



where E can be any proper measure of entanglement of pure states, defined as a monotonically 
increasing function of the entropy of entanglement (i.e. the von Neumann entropy of the 
reduced density matrix of one party). 



If the function E is taken to be exactly the entropy of entanglement, Eq. d53l . then the 
corresponding Gaussian EM is known as Gaussian entanglement of formation [ 106 1 . Any 
Gaussian EM is an entanglement monotone under Gaussian LOCC. The proof given in Sec. 
IV of Ref. 1 106] for the Gaussian entanglement of formation automatically extends to every 
Gaussian EM constructed via the Gaussian convex roof of any proper measure E of pure-state 
entanglement. 

In general, the definition Eq. ( f7TT i involves an optimization over all pure Gaussian states 
with CM cr p smaller than the CM cr of the mixed state whose entanglement one wishes to 
compute. This corresponds to taking the Gaussian convex roof. Despite being a simpler 
optimization problem than that appearing in the definition Eq. (f62b of the true entanglement 
of formation, the Gaussian EMs cannot be expressed in a simple closed form, not even in the 
instance of (nonsymmetric) two-mode Gaussian states. Some recent results on this issue have 
been obtained in 11051 . and will be discussed in section IB~4l in relation with the problem of 
the ordering of quantum states according to different measures of entanglement. 



(70) 



For any mixed Gaussian state g with CM cr, one has [ 

G B (cr)= inf G E (a p ). 




(71) 
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5. Entanglement of two-mode Gaussian states. 

We now discuss the characterization of the prototypical entangled states of CV systems, 
i.e. the two-mode Gaussian states. This includes the explicit determination of the negativities 
and their relationship with global and marginal entropic measures Ml 041 [37l l27l . and the 
evaluation of the Gaussian measures of entanglement 111051 . We will then compare the two 
families of measures according to the ordering that they establish on entangled states. 



5.7. Symplectic parametrization of two-mode Gaussian states. 

To study entanglement and informational properties (like global and marginal entropies) of 
two-mode Gaussian states, we can consider without loss of generality states whose CM er is 
in the Sp(2,n) © Sp(2,B.)-invariant standard form derived in Refs. |[T3ll79ll . Let us recall it 
here for the sake of clarity, 



fa c+ 

a 7 \ a c_ 

7 T P ) = c+ b 

\ c_ b J 



(72) 



For two-mode states, the uncertainty principle Ineq. (fT8l l can be recast as a constraint on 
the 5^(4^) invariants (invariants under global, two-mode symplectic operations) Deter and 
A(cr) = Detc* + Dct/3 + 2Det7 El, 

A(cr) < 1 + Deter . (73) 

The symplectic eigenvalues of a two-mode Gaussian state will be denoted as v_ and v + , 
with V- < v + , with the uncertainty relation (l33l reducing to 



v- > 1 ■ (74) 
A simple expression for the v T can be found in terms of the two Sp( 4 g) invariants II 1001 l28l 



2v\ = A(er) =p y/A 2 (tr) - 4 Deter . (75) 

According to Eq. f72) , two-mode Gaussian states can be classified in terms of their four 
standard form covariances a, b, c+, and c_. It is relevant to provide a reparametrization of 
standard form states in terms of symplectic invariants which admit a direct interpretation for 
generic Gaussian states B104ll37ll27l . Namely, the parameters of Eq. (l72l can be determined 
in terms of the two local symplectic invariants 

Ati = (Dctc*)- 1/2 = l/o, [i 2 = (Dct/3)~ 1/2 = 1/6, (76) 

which are the marginal purities of the reduced single-mode states, and of the two global 
symplectic invariants 

/i = (Det er)- 1 / 2 = [(ab-c^ab-c^)]- 1 / 2 , A = a 2 +fo 2 +2c+c_ , (77) 

which are, respectively, the global purity Eq. d46*l > and the seralian Eq. d32l . Eqs. d76l477] i 
can be inverted to give a physical parametrization of two-mode states in terms of the four 
independent parameters pi, p, and A. This parametrization is particularly useful for the 
evaluation of entanglement ||104|[3"7l . 



Entanglement in continuous variable systems 



25 



5.2. Qualifying and quantifying two-mode entanglement. 

5.2.1. Partial transposition and negativities. The PPT condition for separability, Eq. 
has obviously a very simple form for two-mode Gaussian states. In terms of symplectic 
invariants, partial transposition corresponds to flipping the sign of Det 7, 



(78) 



a 7 \ e _> e T, I a. 7 

7 T fi) ' ^-{l 1 f> 

with Det 7 = —Det 7. For a standard form CM Eq. ( l72l . this simply means c + — > c+, 
c_ — > — c_. Accordingly, the seralian A = Det a + Det (3 + 2 Det 7, Eq. ( 1321 ). is mapped, 
under partial transposition, into 

A = Det a + Det [3 + 2 Det 7 = Det a + Det (3 - 2 Det 7 

= A -4Det 7 = -A + 2/7^ + 2/^2. (79) 

From Eq. (l75l l. the symplectic eigenvalues of the partial transpose <r of a two-mode CM <x are 
promptly determined in terms of symplectic invariants, 



2v% = A t W A 2 - . (80) 



4 

The PPT criterion is then reexpressed by the following inequality 

A<1 + 1/^ 2 , (81) 

equivalent to separability. The state a is separable if and only if i>_ > 1. Accordingly, the 
logarithmic negativity Eq. ( f68b is a decreasing function of z>_ alone, 

Etf = max{0, — log £_} , (82) 

as for the biggest symplectic eigenvalue of the partial transpose one has v + > 1 for all two- 
mode Gaussian states [104, 37]. 

Note that from Eqs. ( I72I73I79I8 lb the following necessary condition for two-mode 
entanglement follows |[T3l . 

cr entangled Det 7 < . (83) 



5.2.2. Entanglement of formation for symmetric states. The optimal convex decomposition 
involved in the definition Eq. d62l of the entanglement of formation [54 1 has been remarkably 
solved in the special instance of two-mode symmetric mixed Gaussian states [i.e. with 
Det a = Det (3 in Eq. d72ll. and turns out to be Gaussian. Namely, the absolute minimum is 
realized within the set of pure two-mode Gaussian states ll97l . yielding 



with 



h{x) 



max [0, h(P-)] 



(1 + xf 
Ax 



log 



(1 



^21 



Ax 



Ax 



log 



(1 



^21 



Ax 



(84) 



(85) 



Such a quantity is, again, a monotonically decreasing function of the smallest symplectic 
eigenvalue i>- of the partial transpose cr of a two-mode symmetric Gaussian CM cr, thus 
providing a quantification of the entanglement of symmetric states equivalent to the one 
provided by the negativities. Lower bounds on the entanglement of formation have been 
derived for nonsymmetric two-mode Gaussian states 01071 . 

As a consequence of the equivalence between negativities and Gaussian measures of 
entanglement in symmetric two-mode Gaussian states, it is tempting to conjecture that there 
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exists a unique quantification of entanglement for all two-mode Gaussian states, embodied 
by the smallest symplectic eigenvalue z>_ of the partially transposed CM, and that the 
different measures simply provide trivial rescalings of the same unique quantification. In 
particular, the ordering induced on the set of entangled Gaussian state is uniquely defined 
for the subset of symmetric two-mode states, and it is independent of the chosen measure of 
entanglement. However, in Sec. 15.41 we will indeed show, within the general framework of 
Gaussian measures of entanglement (see section l4~.5.2l ). that different families of entanglement 
monotones induce, in general, different orderings on the set of nonsymmetric Gaussian states, 
as demonstrated in 1 105 1. 





' a 






' a c_ \ 


lq = ( 


x C+ 




■ 7 P = ( 


v C_ b J 



5.2.3. Gaussian measures of entanglement: geometric framework. The problem of 
evaluating Gaussian measures of entanglement (Gaussian EMs) for a generic two-mode 
Gaussian state has been solved in Ref. 11061 . However, the explicit result contains rather 
cumbersome expressions, involving the solutions of a fourth-order algebraic equation. As 
such, they were judged of no particular insight to be reported explicitly by the authors of 
Ref. Q06). 

We recall here the computation procedure 11051 . For any two-mode Gaussian state 
with CM cr = a s f in standard form Eq. d72l ), a generic Gaussian EM Ge is given by the 
entanglement E of the least entangled pure state with CM er p < cr, see Eq. dTTl i. Denoting by 
7 g (respectively j p ) the 2x2 submatrix obtained from cr by canceling the even (resp. odd) 
rows and columns, we have 

(86) 

All the covariances relative to the "position" operators of the two modes are grouped in j q , 
and analogously for the "momentum" operators in j p . The total CM can then be written as 
a direct sum cr = j q © j p . Similarly, the CM of a generic pure two-mode Gaussian state in 
block-diagonal form (it has been proven that the CM of the optimal pure state has to be with 
all diagonal 2x2 submatrices as well |J06|) can be written as cr p = 7^©7^, where the global 
purity of the state imposes (7^) _1 = j p = T. The pure states involved in the definition of the 
Gaussian EM must thus fulfill the condition 

i; 1 < r < i q . (87) 

This problem admits an interesting geometric formulation II 1 061 - Writing the matrix T 
in the basis constituted by the identity matrix and the three Pauli matrices, 

T= (x +x 3 Xl \ 

y Xl Xq - X 3 ) 

the expansion coefficients (xq, #1,2:3) play the role of coordinates in a three-dimensional 
Minkowski space. In this picture, for example, the rightmost inequality in Eq. ( f8Tb is satisfied 
by matrices T lying on a cone, which is equivalent to the (backwards) light cone of C q in the 
Minkowski space; and similarly for the leftmost inequality. Indeed, one can show that, for the 
optimal pure state cr p opt realizing the minimum in Eq. ( |7T1 i. the two inequalities in Eq. d87b 
have to be simultaneously saturated II 1 0611 . From a geometrical point of view, the optimal T 
has then to be found on the rim of the intersection of the forward and the backward cones of 
jp 1 and 7 9 , respectively. This is an ellipse, and one is left with the task of minimizing the 
entanglement E of cr p — V T -1 (see Eq. d70] i) for T lying on this ellipse[jj| 

|| The geometric picture describing the optimal two-mode state which enters in the determination of the Gaussian 
EMs is introduced in [ 1 06J . A more detailed discussion, including the explicit expression of the Lorentz boost needed 
to move into the plane of the ellipse, can be found in 1 1 08 1 . 
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We recall that any pure two-mode Gaussian state er p is locally equivalent to a two-mode 
squeezed state with squeezing parameter r, described by the CM of Eq. (|2TT >. The following 
statements are then equivalent: (i) E is a monotonically increasing function of the entropy of 
entanglement; (ii) E is a monotonically increasing function of the single-mode determinant 
m 2 = Det a. = Det (3 (see Eq. (f72l); (iii) E is a monotonically decreasing function of the 
local purity /i^ = fi\ = /12 (see Eq. d46b~); (iv) E is a monotonically decreasing function 
of the smallest symplectic eigenvalue v p _ of the partially transposed CM er p ; (v) E is a 
monotonically increasing function of the squeezing parameter r. This chain of equivalences 
is immediately proven by simply recalling that a pure state is completely specified by its 
single-mode marginals, and that for a single-mode Gaussian state there is a unique symplectic 
invariant (the determinant), so that all conceivable entropic quantities are monotonically 
increasing functions of this invariant, as shown in section !?. 3. 21 1 37 1. In particular, statement 
(ii) allows us to minimize directly the single-mode determinant over the ellipse: 

Xi 



= 1 



Detr 



(89) 



with r given by Eq. 

To simplify the calculations, one can move to the plane of the ellipse with a Lorentz boost 
which preserves the relations between all the cones; one can then choose the transformation so 
that the ellipse degenerates into a circle (with fixed radius), and introduce polar coordinates on 
this circle. The calculation of the Gaussian EM for any two-mode Gaussian state is thus finally 
reduced to the minimization of m 2 from Eq. (1891 . at given standard-form covariances of cr, 
as a function of the polar angle 9 on the circle Ml 081 . After some tedious but straightforward 
algebra, one finds 01051 

ml(a,b,c + ,c^) = 

2^1 



1 



c + (ab - <?_) - c_ + cos 9 \J [a - b{ab - c 2 _)\ [b — a(ab — cP_)] 



2(ab-c 2 _) (a 2 + 6 2 + 2c+c_) 



cos 6 [2abc 3 _ + (a 2 



b 2 ) c+c 2 



1 + 2b 2 ) a 2 + b 2 ) c_ - ab (a 2 + b 2 - 2) c+] 



[a — b(ab — c 2 _)\ \b — a(ab — 



(a 2 



[c+(ab 



)+c- 



\ 1 [a - b(ab - (?_)] [b - a{ab - c 2 _)} 



(90) 



where we have assumed c + > |c_| without loss of generality. Therefore, for any entangled 
state, c + > and c_ < 0, see Eq. d83l >. The Gaussian EM (defined in terms of the function 
E on pure states, see Eq. d70l >) of a generic two-mode Gaussian state coincides then with the 
entanglement E computed on the pure state with m 2 — m 2 t , with m 2 t = mine(mg). 



Accordingly, the symplectic eigenvalue i>_ of the partial transpose of the corresponding 
optimal pure-state CM cr p opt , realizing the infimum in Eq. (l7TT i. would read (see Eq. d80l l) 



1. 



-opt — opt) ~ m opt y '"-opt 

As an example, for the Gaussian entanglement of formation II 10611 one has 



G Ef {(t) = h [vl opt (m 2 opt )] 



(91) 



(92) 



with h(x) defined by Eq. d85i 
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Finding the minimum of Eq. d90l > analytically for a generic state is an involved task, but 
numerical optimizations may be quite accurate. 

5.3. Extremal entanglement versus information. 

Here we review the characterization of entanglement of two-mode Gaussian states with 
respect to its relationship with the degrees of information associated with the global state 
of the system, and with the reduced states of each of the two subsystems. 

As extensively discussed in the previous sections, the concepts of entanglement and 
information encoded in a quantum state are closely related. Specifically, for pure states 
bipartite entanglement is equivalent to the lack of information (mixedness) of the reduced 
state of each subsystem. For mixed states, each subsystem has its own level of impurity, 
and moreover the global state is itself characterized by a nonzero mixedness. Each of these 
properties can be interpreted as information on the preparation of the respective (marginal and 
global) states, as clarified in section 13.3.11 Therefore, by properly accessing these degrees 
of information it should be possible to deduce, at least to some extent, the status of the 
correlations between the subsystems. 

Indeed, the original studies revealed that, at fixed global and marginal degrees of purity 
(or of generalized entropies), the negativities of arbitrary (mixed) two-mode Gaussian states 
are analytically constrained by rigorous upper and lower bounds Ill04ll37ll27l . This follows 
by reparametrizing, as already anticipated, the standard form CM Eq. (l72b in terms of 
the invariants /j,2, /i, A and by observing that, at fixed purities, the negativities are 
monotonically decreasing function of A. Further constraints imposed on A by the uncertainty 
principle and by the existence condition of the radicals involved in the reparametrization, 



Gaussian states at fixed global and local purities. They are known, respectively, as "GMEMS" 
(saturating the leftmost inequality in Eq. d93l). alias nonsymmetric thermal squeezed states, 
and "GLEMS" (saturating the rightmost inequality in Eq. d93l). alias mixed states of partial 
minimum uncertainty II 1041 [37). Nonsymmetric thermal squeezed states have also been 
proven to be maximally entangled Gaussian mixed states at fixed global purity and mean 
energy |fT09l . 

Summarizing, the entanglement, quantified by the negativities, of two-mode ( mixed) 
Gaussian states is strictly bound from above and from below by the amounts of global 
and marginal purities, with only one remaining degree of freedom related to the symplectic 
invariant A. 

The existence of GMEMS and GLEMS has two consequences. First, it allows for a 
classification of the properties of separability of all two-mode Gaussian states according 
to their degree of global and marginal purities, as summarized in Table @] Namely, from 
the separability properties of the extremally entangled states, necessary and/or sufficient 
conditions for entanglement — which constitute the strongest entropic criteria for separability 
1 1 10 1 to date in the case of Gaussian states — are straightforwardly derived, which allow one 
to decide almost unambiguosly if a given two-mode Gaussian state is entangled or not based 
on its degree of global and local purities. There is only a narrow region where, for given 
purities, both separable and entangled states can coexist, as pictorially shown in Fig.Q] 

The second consequence is of a quantitative nature. The quantitative analysis of the 
maximal and the minimal entanglement allowed for a Gaussian state with given purities, 




(93) 



immediately lead to the definition of extremally - maximally and minimally - entangled 
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Figure 1. Summary of entanglement properties of two-mode (nonsymmetric) Gaussian states 
in the space of marginal purities and global purity fi (we plot the normalized ratio 

Ii/filfi2 to gain a better graphical clarity). In this space, separable states (red zone) and 
entangled states (green to magenta zone, according to the average entanglement) are well 
separated except for a narrow region of coexistence (depicted in yellow). The mathematical 
relations defining the boundaries between the three regions are collected in Table|4] The three- 
dimensional envelope is cut at z = 3.5. 

reveals that they are very close to each other, their difference narrowing exponentially 
with increasing entanglement. One can then define the average logarithmic negativity 
(mean value of the entanglements of the GMEMS and the GLEMS corresponding to a 
given triplet of purities) as a reliable estimator of bipartite entanglement in two-mode 
Gaussian states and its accurate quantification by knowledge of the global and marginal 
purities alone II 1 0411 . The purities are nonlinear functionals of the state, but, assuming 
some prior knowledge about it (essentially, its Gaussian character), be measured through 
direct methods, in particular by means of single-photon detection schemes Mil 11 (of which 
preliminary experimental verifications are available II 121 ) or by interferometric quantum- 
network architectures 111 131 II 141 II 151 . Very recently, a scheme to measure locally all 
symplectic invariants (and hence the entanglement) of two-mode Gaussian states has been 
proposed [116|. Notice that no complete homodyne reconstruction |117| of the covariance 



Degrees of purity 


Entanglement properties 




unphysical region 


UlU 2 < u < J 11 ^ 2 


separable states 


MiA»2 ^ it ^ ' tlM2 

M1+P2-P1M2 ^ '~ - y/nl + fi-ptnl 


coexistence region 


M1M2 ^ u ^ ** 1 ' <a 
y/lil+Ma-Hil4 ' ~ M1M2 + IM1-M2I 


entangled states 


U > M ^ M2 

^1^2 + 1/^1-^2 


unphysical region 



Table 4. Classification of two-mode Gaussian states and of their properties of separability 
according to their degrees of global purity fi and of marginal purities /11 and fi2 II 0411371 . 
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matrix is needed in all these schemes. 

The average estimate of the logarithmic negativity becomes indeed an exact 
quantification in the two important instances of GMEMS (nonsymmetric thermal squeezed 
states) and GLEMS (mixed states of partial minimum uncertainty), whose logarithmic 
negativity is completely determined as a function of the three purities alone II 1 041 [371 - 

5.4. Ordering Gaussian states with measures of entanglement. 

The role of GMEMS and GLEMS in the characterization of entanglement of two-mode 
Gaussian states is not confined to the choice of the negativities. In fact, the Gaussian EMs for 
these two families of states can be obtained in closed form, as the optimization involved in 
Eq. d90l ) admits a rather simple analytical solutions in this two instances, as we have shown in 
Ref. [105]. One might therefore raise the question whether such extremally entangled states 
(with respect to the negativities) conserve their role once the entanglement is measured by 
the Gaussian EMs. The answer, somehow surprisingly, is a negative one: actually, there is 
a large region in the space of purities, where GMEMS are strictly less entangled than the 
corresponding GMEMS, when their entanglement is measured according to the Gaussian 
EMs. 

The Gaussian EMs and the negativities are thus not equivalent for the quantification of 
entanglement in mixed, nonsymmetric two-mode Gaussian states B105II . The interpretation 
of this result needs to be discussed with some care. On the one hand, one could think 
that the ordering induced by the negativities is a natural one, due to the fact that such 
measures of entanglement are directly inspired by the necessary and sufficient PPT criterion 
for separability. Thus, one would expect that the ordering induced by the negativities should 
be preserved by any bona fide measure of entanglement, especially if one considers that 
the extremal states, GLEMS and GMEMS, have a clear physical interpretation. Therefore, 
as the Gaussian entanglement of formation is an upper bound to the true entanglement of 
formation, one could be tempted to take this result as an evidence that the latter is globally 
minimized on non Gaussian decomposition, at least for GLEMS. However, this is only a 
qualitative/speculative argument: proving or disproving that the Gaussian entanglement of 
formation is the true one for any two-mode Gaussian state remains an open question [56]. 

On the other hand, one could take the simplest discrete-variable instance, constituted by a 
two-qubit system, as a test-case for comparison. There, although for pure states the negativity 
coincides with the concurrence, an entanglement monotone equivalent to the entanglement 
of formation for all states of two qubits 111 1 81 . the two measures cease to be equivalent for 
mixed states, and the orderings they induce on the set of entangled states can be different 
111 191 11201 . This analogy seems to support the stand that, in the arena of mixed states, the 
definition of a unique measure of entanglement cannot really be pursued, due to the different 
operational meanings and physical processes (in the cases when it has been possible to identify 
them) that are associated to each definition: one could think, for instance, of the operational 
difference existing between the definitions of distillable entanglement and entanglement cost 
lf54ll . In other words, from this point of view, each inequivalent measure of entanglement 
introduced for mixed states should capture physically distinct aspects and forms of quantum 
correlations existing in these states. Gaussian EMs might then still be considered as proper 
measures of CV entanglement, adapted to a different context than negativities. This view 
seems particularly appropriate when constructing Gaussian EMs to investigate entanglement 
sharing in multipartite Gaussian states 111211 11221. as discussed in section|7] 

The inequivalence between the two families of CV entanglement measures is somehow 
tempered. Namely, we have rigorously proven that, at fixed negativities, the Gaussian 
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Figure 2. Comparison between the Gaussian entanglement of formation Ge f and the 
logarithmic negativity E^f for two-mode Gaussian states. Symmetric states accommodate 
on the lower boundary (solid line). States of maximal negativities at fixed (infinite) average 
local mixedness, lie on the dashed line. All GMEMS and GLEMS lie below the dashed line. 
The latter is conjectured, with strong numerical support, to be the upper boundary for the 
Gaussian entanglement of formation of all two-mode Gaussian states, at fixed negativity. Cfr. 
Ref. 1 105 1 for the complete discussion and the detailed mathematical proofs. 

measures of entanglement are bounded from below (the states which saturate this bound are 
simply symmetric two-mode states); moreover, we provided some strong evidence suggesting 
that they are as well bounded from above B105II . A direct comparison between the logarithmic 
negativity and the Gaussian entanglement of formation for all two-mode Gaussian states is 
shown in Fig. [2] 

6. Multimode bipartite entanglement: localization and scaling. 

We shall now discuss the properties of entanglement in multimode Gaussian states endowed 
with particular symmetry constraints under permutations of the modes ll46l |45ll . The 
usefulness of these states arise in contexts like quantum error correction 11231 , where some 
redundancy is required for a fault-tolerant encoding of information. Bisymmetric and, as a 
special case, fully symmetric Gaussian states have been introduced in section [J.4.2l The study 
of the symplectic spectra of (M + 7V)-mode Gaussian states reveals that, with respect to the 
bipartition across which they exhibit the local permutation invariance (any bipartition is valid 
for fully symmetric states), local symplectic diagonalizations of the M -mode and the TV-mode 
blocks result in a complete reduction of the multimode state to an equivalent two-mode state, 
tensored with M + N — 2 uncorrected thermal single-mode states. The equivalent two- 
mode state possesses all the information of the original bisymmetric multimode state for what 
concerns entropy and entanglement. As a consequence, the validity of the PPT criterion as a 
necessary and sufficient condition for separability has been extended to bisymmetric Gaussian 
states in section l4.4.1l 

6.1. Unitarily localizable entanglement of Gaussian states. 

Here, equipped with the tools introduced for the analysis of two-mode entanglement in 
Gaussian states, we discuss multimode entanglement in symmetric and bisymmetric Gaussian 
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states. In particular, we will investigate how the block entanglement scales with the number of 
modes at fixed squeezing. The form of the scaling hints at the presence of genuine multipartite 
entanglement that progressively arises among all the modes as their total number increases. 

The central observation of the present section is contained in the following result 
11461 l45ll . straightforwardly deducible from the discussions in section 13.4.21 and section 
14.4.11 The bipartite entanglement of bisymmetric (M + N)-mode Gaussian states under 
M X N partitions is "unitarily localizable", namely, through local unitary (reversible) 
operations, it can be completely concentrated onto a single pair of modes, each of them 
belonging respectively to the M-mode and to the N-mode blocks. Hence the multimode 
block entanglement (i.e. the bipartite entanglement between two blocks of modes) of 
bisymmetric (generally mixed) Gaussian states can be determined as an equivalent two-mode 
entanglement. The entanglement will be quantified by the logarithmic negativity in the general 
instance because the PPT criterion holds, but we will also show some explicit nontrivial cases 
in which the entanglement of formation, Eq. d62l . between M -mode and iV-mode parties can 
be exactly computed. 

We remark that the notion of "unitarily localizable entanglement" is different from that 
introduced by Verstraete, Popp, and Cirac for spin systems [124]. There, it was defined 
as the maximal entanglement that can be concentrated on two chosen spins through local 
measurements on all the other spins. Here, the local operations that concentrate all the 
multimode entanglement on two modes are unitary and involve the two chosen modes as 
well, as parts of the respective blocks. The localizable entanglement in the sense of [ 124| can 
be computed as well for (mixed) symmetric Gaussian states of an arbitrary number of modes, 
and is in direct quantitative connection with the optimal fidelity of multiparty teleportation 
networks 111251 . 

It is important to observe that the unitarily localizable entanglement (when computable) 
is always stronger than the localizable entanglement in the sense of 01241 . In fact, if we 
consider a generic bisymmetric multimode state of a M x N bipartition, with each of the 
two target modes owned respectively by one of the two parties (blocks), then the ensemble of 
optimal local measurements on the remaining ("assisting") M + N — 2 modes belongs to the 
set of LOCC with respect to the considered bipartition. By definition the entanglement cannot 
increase under LOCC, which implies that the localized entanglement (in the sense of II124II ) 
is always less or equal than the original M x N block entanglement. On the contrary, all of 
the same M x N original bipartite entanglement can be unitarily localized onto the two target 
modes. 

This is a key point, as such local unitary transformations are reversible by definition. 
Therefore, by only using passive and active linear optics elements such as beam splitters, 
phase shifters and squeezers, one can in principle implement a reversible machine 
(entanglement switch) that, from mixed, bisymmetric multimode states with strong quantum 
correlations between all the modes (and consequently between the A/-mode and the N- 
mode partial blocks) but weak couplewise entanglement, is able to extract a highly pure, 
highly entangled two-mode state (with no entanglement loss, as all the M x N entanglement 
can be localized). If needed, the same machine would be able, starting from a two-mode 
squeezed state and a collection of uncorrelated thermal or squeezed states, to distribute 
the two-mode entanglement between all modes, converting the two-mode into multimode, 
multipartite quantum correlations, again with no loss of entanglement. The bipartite or 
multipartite entanglement can then be used on demand, the first for instance in a CV quantum 
teleportation protocol 11261 . the latter to secure quantum key distribution or to perform 
multimode entanglement swapping 11271 . We remark, once more, that such an entanglement 
switch is endowed with maximum (100%) efficiency, as no entanglement is lost in the 
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conversions. 



6.2. Quantification and scaling of entanglement in fully symmetric states. 

Here we will review the explicit evaluation of the bipartite block entanglement for some 
instances of multimode Gaussian states. We will discuss its scaling behavior as a function 
of the number of modes and explore in some detail the localizability of the multimode 
entanglement. We focus our attention on fully symmetric L-mode Gaussian states (the number 
of modes is denoted by L in general to avoid confusion), endowed with complete permutation 
invariance under mode exchange, and described by a 2L x 2L CM ctrl given by Eq. ( Bil l. 
These states are trivially bisymmetric under any bipartition of the modes, so that their block 
entanglement is always localizable by means of local symplectic operations. Let us recall that 
concerning the covariances in normal forms of fully symmetric states (see section[3A2]l, pure 
L-mode states are characterized by 

_ (L-2)(fr 2 -l) + A /(62-l)[L((fe2-l)L + 4) = 4l 
Zl 26(L-1) 

(94) 

(L - 2)(b 2 - 1) - ^(fr 2 - 1) [L ((b 2 - 1) L + 47^4] 
Z2 26(L-1) 

Pure, fully symmetric Gaussian states are generated as the outputs of the application of a 
sequence of L — 1 beam splitters to L single-mode squeezed inputs. 1112811901 . The CM cr p ^ L 
of this class of pure states, for a given number of modes, depends only on the parameter 
b = l/fip > 1, which is an increasing function of the single-mode squeezing needed 
to prepare the state. Correlations between the modes are induced according to the above 
expression for the covariances Their multipartite entanglement sharing has been studied 
in Ref. 11211 (see next section), and their use in teleportation networks have been investigated 
inRefs. |[T28l[T25l . 

In general, one can compute the entanglement between a block of K modes and the 
remaining L — K modes for pure states (in this case the block entanglement is simply 
equivalent to the von Neumann entropy of each of the reduced blocks) and for mixed, fully 
symmetric states under any bipartition of the modes. 



6.2.7. 1 x iV entanglement. Based on Ref. Il46l . we begin by assigning a single mode to 
subsystem Sa, and an arbitrary number N of modes to subsystem Sb, forming a CV system 
globally prepared in a fully symmetric (1 + iV)-mode Gaussian state of modes. 

We consider pure, fully symmetric states with CM er^ 1+JV , obtained by inserting Eq. d94l i 
with L = (1 + N). Exploiting our previous analysis, we can compute the entanglement 
between a single mode with reduced CM cr 13 and any i-C-mode partition of the remaining 

@\8 K 

modes (1 < K < N), by determining the equivalent two-mode CM creq ■ We remark that, 
for every K, the 1 x K entanglement is always equivalent to a 1 x 1 entanglement, so that 
the quantum correlations between the different partitions of cr can be directly compared to 
each other: it is thus possible to establish a multimode entanglement hierarchy without any 
problem of ordering. 

0\0 K 

The 1 x K entanglement quantified by the logarithmic negativity E^j- is determined 
by the smallest symplectic eigenvalue C/^' N ^ of the partially transposed CM a^J 3 . For any 
nonzero squeezing (i.e. b > 1) one has that f/^' N > < l, meaning that the state exhibits 
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Figure 3. Entanglement hierarchy for (1 + A r )-mode fully symmetric pure Gaussian states 
(AT = 9). 




Figure 4. Scaling as a function of N of the lxl entanglement (green bars) and of the 
1 X (N — 1) entanglement (red bars) for a (1 + A f )-mode pure fully symmetric Gaussian 
state, at fixed squeezing (b = 1.5). 



genuine multipartite entanglement: each mode is entangled with any other A'-mode block, as 
first remarked in Ref. ]90ll . Further, the genuine multipartite nature of the entanglement can 
be precisely quantified by observing that 

E M > E M , 

as shown in Fig. [3] 

The lxl entanglement between two modes is weaker than the 1x2 one between a mode 
and other two modes, which is in turn weaker than the 1 x K one, and so on with increasing 
K in this typical cascade structure. From an operational point of view, a signature of genuine 
multipartite entanglement is revealed by the fact that performing e.g. a local measurement on 
a single mode will affect all the other N modes. This means that the quantum correlations 
contained in the state with CM <t^ 1+n can be fully recovered only when considering the 1 x N 
partition. 

In particular, the pure-state 1 x N logarithmic negativity is, as expected, independent of 
N, being a simple monotonic function of the entropy of entanglement Ey, Eq. d53l l (defined 
as the von Neumann entropy of the reduced single-mode state with CM a$). It is worth noting 
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that, in the limit of infinite squeezing (b — > oo), only the 1 x N entanglement diverges while 
all the other 1 x K quantum correlations remain finite (see Fig. [3j. Namely, 



uP\P K ( p \ b^oo 1 . 



1-4K 



(95) 



_N(K + 1) -K(K~3) 

which cannot exceed log ~ 0.8 for any N and for any K < N. 

At fixed squeezing (i.e. fixed local properties, b = 1/np), the scaling with N of the 
1 x (N — 1) entanglement compared to the lxl entanglement is shown in Fig.|4](we recall 
that the 1 x N entanglement is independent on N). Notice how, with increasing number of 
modes, the multimode entanglement increases to the detriment of the two-mode one. The 
latter is indeed being distributed among all the modes: this feature will be properly quantified 
within the framework of CV entanglement sharing in the next section II121I . 

We remark that such a scaling feature occurs both in fully symmetric and bisymmetric 
states (think, for instance, to a single-mode squeezed state coupled with a TV-mode symmetric 
thermal squeezed state), pure or mixed. The simplest example of a mixed state in which our 
analysis reveals the presence of genuine multipartite entanglement is obtained from c^i+w 
by tracing out some of the modes. Fig. [4] can then also be seen as a demonstration of the 
scaling in such a iV-mode mixed state, where the 1 x (N — 1) entanglement is the strongest 
one. Thus, with increasing N, the global mixedness can limit but not destroy the distribution 
of entanglement in multiparty form among all the modes. 



6.2.2. M x N entanglement. Based on the results of Ref. P31 , one can consider a generic 
27V-mode fully symmetric mixed state with CM <t p ^, see Eq. (IBTT l. obtained from a pure 
fully symmetric (2N + Q)-mode state by tracing out Q modes. 

For any Q, for any dimension K of the block (K < N), and for any non zero 
squeezing (i.e. for b > 1) one has that vk < 1, meaning that the state exhibits genuine 
multipartite entanglement, generalizing the lxAf case described before: each K-mode party 
is entangled with the remaining (27V — if) -mode block. The genuine multipartite nature of 

the entanglement can be determined by observing that, again, E\f is an increasing 

function of the integer K < N, as shown in Fig. [5] The multimode entanglement of mixed 
states remains finite also in the limit of infinite squeezing, while the multimode entanglement 
of pure states diverges with respect to any bipartition, as shown in Fig. [5] 

In fully symmetric Gaussian states, the block entanglement is unitarily localizable with 
respect to any K x (2N — K) bipartition. Since in this instance all the entanglement can 
be concentrated on a single pair of modes, after the partition has been decided, no strategy 
could grant a better yield than the local symplectic operations bringing the reduced CMs in 
Williamson form (because of the monotonicity of the entanglement under general LOCC). 
However, the amount of block entanglement, which is the amount of concentrated two-mode 
entanglement after unitary localization has taken place, actually depends on the choice of a 
particular K x (2N — K) bipartition, giving rise to a hierarchy of localizable entanglements. 

Let us imagine that for a given Gaussian multimode state (say, for simplicity, a fully 
symmetric state) its entanglement is meant to serve as a resource for a given protocol. Let us 
next suppose that the protocol is optimally implemented if the entanglement is concentrated 
between only two modes of the global systems, as it is the case, e.g., in a CV teleportation 
protocol between two single-mode parties [126]. Which choice of the bipartition between 
the modes allows for the best entanglement concentration by a succession of local unitary 
operations? In this framework, it turns out that assigning K = 1 mode at one party and all the 
remaining modes to the other, as discussed in section |6". 2. II constitutes the worst localization 
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Figure 5. Hierarchy of block entanglements of fully symmetric 27V-mode Gaussian states 
of K X (27V — K) bipartitions (27V = 10) as a function of the single-mode squeezing b. 
The block entanglements are depicted both for pure states (solid lines) and for mixed states 
obtained from fully symmetric (27V + 4)-mode pure Gaussian states by tracing out 4 modes 
(dashed lines). 



strategy B31 . Conversely, for an even number of modes the best option for localization is an 
equal K = N splitting of the 2N modes between the two parties. The logarithmic negativity 

N \B N 

E\f , concentrated into two modes by local operations, represents the optimal localizable 
entanglement (OLE) of the state (Tqin, where "optimal" refers to the choice of the bipartition. 

B N + 1 \B N 

Clearly, the OLE of a state with 2TV + 1 modes is given by E^ . These results may be 
applied to arbitrary, pure or mixed, fully symmetric Gaussian states. 

We now turn to the study of the scaling behavior of the OLE of 27V-mode Gaussian states 
when the number of modes is increased, to understand how the number of local cooperating 
parties can improve the maximal entanglement that can be shared between two parties. For 
generic (mixed) fully symmetric 2iV-mode states of N x N bipartitions, the OLE can be 
quantified also by the entanglement of formation Ep, Eq. (|84| >. as the equivalent two-mode 
state is symmetric ll45l . It is then useful to compare, as a function of TV, the lxl entanglement 
of formation between a pair of modes (all pairs are equivalent due to the global symmetry of 
the state) before the localization, and the N x TV entanglement of formation, which is equal 
to the optimal entanglement concentrated in a specific pair of modes after performing the 
local unitary operations. The results of this study are shown in Fig. [6] The two quantities 
are plotted at fixed squeezing 6 as a function of N both for a pure 2V-mode state with CM 
cr^ 2JV and a mixed 2n-mode state with CM c^t . As the number of modes increases, any 
pair of single modes becomes steadily less entangled, but the total multimode entanglement 
of the state grows and, as a consequence, the OLE increases with V. In the limit V — > oo, 
the N x N entanglement diverges while the lxl one vanishes. This exactly holds both for 
pure and mixed states, although the global degree of mixedness produces the typical behavior 
that tends to reduce the total entanglement of the state. 

6.3. Towards multipartite entanglement. 

We have shown that bisymmetric multimode Gaussian states (pure or mixed) can be reduced, 
by local symplectic operations, to the tensor product of a correlated two-mode Gaussian 
state and of uncorrelated thermal states (the latter being obviously irrelevant as far as the 
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Figure 6. Scaling, with half the number of modes, of the entanglement of formation in two 
families of fully symmetric 27V-mode Gaussian states. Plot (a) depicts pure states, while 
mixed states (b) are obtained from (2N + 4)-mode pure states by tracing out 4 modes. For 
each class of states, two sets of data are plotted, one referring to N X N entanglement (yellow 
bars), and the other to 1 X 1 entanglement (blue bars). Notice how the N X N entanglement, 
equal to the optimal localizable entanglement (OLE) and estimator of genuine multipartite 
quantum correlations among all the 2N modes, increases at the detriment of the bipartite 
lxl entanglement between any pair of modes. The single-mode squeezing parameter is fixed 
at b = 1.5. 

correlation properties of the multimode Gaussian state are concerned). As a consequence, all 
the entanglement of bisymmetric multimode Gaussian states of arbitrary M x N bipartitions 
is unitarily localizable in a single (arbitrary) pair of modes shared by the two parties. Such 
a useful reduction to two-mode Gaussian states is somehow similar to the one holding for 
states with fully degenerate symplectic spectra [40, 41 J, encompassing the relevant instance 
of pure states, for which all the symplectic eigenvalues are equal to 1. The present result 
allows to extend the PPT criterion as a necessary and sufficient condition for separability 
for all bisymmetric multimode Gaussian states of arbitrary M x N bipartitions (as shown in 
section |4~.4.1| i, and to quantify their entanglement |46l l45l . 

Notice that, in the general bisymmetric instance addressed in this section, the possibility 
of performing a two-mode reduction is crucially partition-dependent. However, as we have 
explicitly shown, in the case of fully symmetric states all the possible bipartitions can be 
analyzed and compared, yielding remarkable insight into the structure of the multimode block 
entanglement of Gaussian states. This leads finally to the determination of the maximum, or 
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optimal localizable entanglement that can be concentrated on a single pair of modes. 

The multipartite entanglement in the considered class of multimode Gaussian states can 
be produced, detected [128, 90 1, and, by virtue of the present analysis, reversibly localized 
by all-optical means. Moreover, the multipartite entanglement allows for a reliable {i.e. with 
fidelity .F > Td, where T c i = 1/2 is the classical threshold 1 129l [T30lD quantum teleportation 
between any two parties with the assistance of the remaining others B128L The connection 
between entanglement in the symmetric Gaussian resource states and optimal teleportation- 
network fidelity has been clarified in 111251 . as briefly recalled in section [9721 

The present section is intended as a bridge between the two main parts of this 
review, the one dealing with bipartite entanglement and the one dealing with multipartite 
entanglement. We have characterized entanglement in multimode Gaussian states by reducing 
it to a two-mode problem. By comparing the equivalent two-mode entanglements in the 
different bipartitions one is able to detect unambiguously the presence of genuine multipartite 
entanglement. We now analyze in more detail the sharing phenomenon responsible for 
the distribution of entanglement from a bipartite, two-mode form, to a genuine multipartite 
manifestation in iV-mode Gaussian states, both with and without symmetry constraints. 

7. Distributed entanglement and monogamy inequality for all Gaussian states. 

We will now review some recent work on multipartite entanglement sharing in CV systems, 
that has lead to the definition of a mathematically and physically bona fide measure of genuine 
tripartite entanglement for arbitrary three-mode Gaussian states 1121111311 , the proof of the 
monogamy inequality on distributed entanglement for all Gaussian states [122|, and the 
demonstration of the promiscuous sharing structure of multipartite entanglement in Gaussian 
states 01211 . which can be unlimited in states of more than three modes [ 132|. 

We begin by introducing some new entanglement monotones, the contangle, the 
Gaussian contangle, and the Gaussian tangle II121II1221 . apt to quantify distributed Gaussian 
entanglement, thus generalizing to the CV setting the tangle [133] defined for systems of 
qubits. 

Motivated by the analysis of the block entanglement hierarchy and its scaling structure 
in fully symmetric Gaussian states (see previous section) we will review the central result that 
CV entanglement, once properly quantified, is monogamous for all Gaussian states, in the 
sense that it obeys a proper generalization to CV systems of the Coffman-Kundu-Wootters 
monogamy inequality 11221 . In the next section we will review the results concerning the 
simplest instance of tripartite CV entangled states, the three-mode Gaussian states 11311 . 
For these states, thanks to the monogamy inequality, it is possible to contruct a measure 
of genuine tripartite entanglement, the residual contangle, that has been shown to be a 
full entanglement monotone under Gaussian LOCC 111211 . Equipped with such a powerful 
tool to quantify tripartite entanglement, we will proceed to review the entanglement sharing 
structure in three-mode Gaussian states, and the related, peculiar property of promiscuity for 
CV entanglement. This property essentially consists in the fact that bipartite and multipartite 
entanglement in multimode Gaussian states can be enhanced and simultaneously maximized 
without violating the monogamy inequality on entanglement sharing, and can even grow 
unboundedly in Gaussian states of more than three modes 111321 . 

7.1. The need for a new continuous-variable entanglement monotone. 

Our primary aim, as in Ref. 111211 . is to analyze the distribution of entanglement between 
different (partitions of) modes in Gaussian states of CV systems. 
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In Ref. 111331 Coffman, Kundu and Wootters (CKW) proved for system of three qubits, 
and conjectured for N qubits (this conjecture has now been proven by Osborne and Verstraete 
[ 134 1), that the bipartite entanglement E (properly quantified) between, say, qubit A and the 
remaining two-qubits partition (BC) is never smaller than the sum of the A|B and A|C bipartite 
entanglements in the reduced states: 

E A\{BC) > E A\B +E A\C_ (96) 

This statement quantifies the so-called monogamy of quantum entanglement 11351 . in 
opposition to the classical correlations, which are not constrained and can be freely shared 
Ifl36l . 

One would expect a similar inequality to hold for three-mode Gaussian states, namely 

£«I(J*) _ J5«|3 _ £«|fc > , (97) 

where E is a proper measure of bipartite CV entanglement and the indexes {i, j, k} label the 
three modes. However, the demonstration of such a property is plagued by subtle difficulties. 

Let us consider the simplest conceivable instance of a pure three-mode Gaussian state 
completely invariant under mode permutations. These pure Gaussian states are named fully 
symmetric (see section 13.4.2b . and their standard form CM (obtained by inserting Eq. d9"4i > 
with L = 3 into Eq. f5i\ ) is only parametrized by the local mixedness b = (l//tts) > 1> an 
increasing function of the single-mode squeezing r; oc , with b — > 1 + when r; oc — * + . For 
these states, it is not difficult to show that the inequality (|97l i can be violated for small values 
of the local squeezing factor, using either the logarithmic negativity £V or the entanglement 
of formation Ep (which is computable in this case via Eq. (1841 1. because the two-mode 
reduced mixed states of a pure symmetric three-mode Gaussian states are again symmetric) to 
quantify the bipartite entanglement. This fact implies that none of these two measures is the 
proper candidate for approaching the task of quantifying entanglement sharing in CV systems. 
This situation is reminiscent of the case of qubit systems, for which the CKW inequality 
holds using the tangle r, defined as the square of the concurrence [118], but can fail if one 
chooses equivalent measures of bipartite entanglement such as the concurrence itself or the 
entanglement of formation 11331 . 

It is then necessary to define a proper measure of CV entanglement that specifically 
quantifies entanglement sharing according to a monogamy inequality of the form (|97l i 11211 . 
A first important hint toward this goal comes by observing that, when dealing with 1 x N 
partitions of fully symmetric multimode pure Gaussian states together with their lxl 
reduced partitions, the desired measure should be a monotonically decreasing function / of 
the smallest symplectic eigenvalue i>- of the corresponding partially transposed CM <x. This 
requirement stems from the fact that i>- is the only eigenvalue that can be smaller than 1, 
as shown in section 14.4.11 violating the PPT criterion with respect to the selected bipartition. 
Moreover, for a pure symmetric three-mode Gaussian state, it is necessary to require that 
the bipartite entanglements E l ^ k ^ and E % ^ — E^ k be respectively functions f{n^ k ^) and 

f{v % -) of the associated smallest symplectic eigenvalues z)^ and , in such a way that 
they become infinitesimal of the same order in the limit of vanishing local squeezing, together 
with their first derivatives: 

K - .,. = y - , . -> 1 for b -► 1+ , (98) 
2f{v lb ) 2/'(^) 

where the prime denotes differentiation with respect to the single-mode mixedness b. The 
violation of the sharing inequality d97l ) exhibited by the logarithmic negativity can be in fact 
traced back to the divergence of its first derivative in the limit of vanishing squeezing. The 
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above condition formalizes the physical requirement that in a symmetric state the quantum 
correlations should grow smoothly and be distributed uniformly among all the three modes. 
One can then see that the unknown function / exhibiting the desired property is simply the 
squared logarithmic negativity^ 

/(*_) = [-log^] 2 . (99) 

We remind again that for (fully symmetric) (1 + iV)-mode pure Gaussian states, the partially 
transposed CM with respect to any 1 x JV bipartition, or with respect to any reduced lxl 
bipartition, has only one symplectic eigenvalue that can drop below 1; hence the simple form 
of the logarithmic negativity (and, equivalently, of its square) in Eq. (|99l . 



7.2. Squared negativities as continuous-variable tangles. 

Equipped with this finding, one can give a formal definition of a bipartite entanglement 
monotone E T that, as we will soon show, can be regarded as a CV analogue of the tangle. 
Note that the context here is completely general and we are not assuming that we are dealing 
with Gaussian states only. For a generic pure state of a (1 + iV)-mode CV system, we 
define the square of the logarithmic negativity (the latter defined by Eq. d66l>): 

■®t(V0 = i°g 2 II fill l > Q = \1>M- (ioo) 

This is a proper measure of bipartite entanglement, being a convex, increasing function of 
the logarithmic negativity Ej\f, which is equivalent to the entropy of entanglement Eq. d53l 
for arbitrary pure states in any dimension. Def. dlOOb is naturally extended to generic mixed 
states p of (N + l)-mode CV systems through the convex-roof formalism. Namely, we can 
introduce the quantity 

E T {p) = inf V ft £ T (i) , (101) 

where the infimum is taken over all convex decompositions of p in terms of pure states {IV'i)}* 
and if the index i is continuous, the sum in Eq. ( 11011 ) is replaced by an integral, and the 
probabilities {pi} by a probability distribution Let us now recall that, for two qubits, 

the tangle can be defined as the convex roof of the squared negativity Ml 371 (the latter being 
equal to the concurrence Ml 1811 for pure two-qubit states 11201 11371 . Here, Eq. ( 11011 ) states 
that the convex roof of the squared logarithmic negativity properly defines the continuous- 
variable tangle, or, in short, the contangle E T (p), in which the logarithm takes into account 
for the infinite dimensionality of the underlying Hilbert space. 

On the other hand, by recalling the equivalence of negativity and concurrence for pure 
states of qubits, the tangle itself can be defined for CV systems as the convex-roof extension 
of the squared negativity. Let us recall that the negativity JV, Eq. d63l of a quantum state q is 
a convex function of the logarithmic negativity Etf, Eq. d66T ). 

JV(g) = eXp[ ^ (g) = 1] . (102) 

Notice that an infinite number of functions satisfying Eq. d98t can be obtained by expanding f{p— ) around z>_ = 1 
at any even order. However, they are all monotonic convex functions of /. If the inequality )97t holds for /, it will 
hold as well for any monotonically increasing, convex function of /, such as the logarithmic negativity raised to any 
even power k > 2, but not for ft = 1. We will exploit this "gauge freedom" in the following, to define an equivalent 
entanglement monotone in terms of squared negativity 1 122|. 
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7.2.7. Gaussian contangle and Gaussian tangle. From now on, we will restrict our attention 
to Gaussian states. 

Gaussian contangle. — For any pure multimode Gaussian state with CM cr p , of 
N + 1 modes assigned in a generic lxJV repartition, explicit evaluation gives immediately 
that E T (ip) = E T (cr p ) takes the form 

E T (o*) = log 2 (l/vi - y/l/fi ~ X ) - ( 103 ) 

where [i\ = 1/VDct ct\ is the local purity of the reduced state of mode 1 with CM <j\ . 

For any multimode, mixed Gaussian states with CM cr, we will then denote the contangle 
by E T (cr), in analogy with the notation used for the contangle E T (cr p ) of pure Gaussian 
states in Eq. (11031 . Any multimode mixed Gaussian state with CM cr, admits at least 
one decomposition in terms of pure Gaussian states er p only. The infimum of the average 
contangle, taken over all pure Gaussian state decompositions, defines then the Gaussian 
contangle G T : 

G T (a) = inf / n(dcr p )E T (cr p ) . (104) 

{ir(dirP),CTP} J 

It follows from the convex roof construction that the Gaussian contangle G T (cr) is an upper 
bound to the true contangle E T (cr) (as the latter can be in principle minimized over a non 
Gaussian decomposition), 

E T (cr) < G T (cr) , (105) 

It can be shown that G T (cr) is a bipartite entanglement monotone under Gaussian LOCC: 
in fact, the Gaussian contangle belongs to the general family of Gaussian entanglement 
measures, whose properties as studied in Ref. [105] have been presented in section 
14.5.21 Therefore, for Gaussian states, the Gaussian contangle, similarly to the Gaussian 
entanglement of formation 11061 , takes the simple form 

G r (cr) = inf E T (cr p ) , (106) 

crP<cr 

where the infimum runs over all pure Gaussian states with CM <j p < cr. Let us remark that, if 
cr denotes a mixed symmetric two-mode Gaussian state, then the Gaussian decomposition is 
the optimal one [97] (see section l5~.2.2b . and the optimal pure-state CM <j p minimizing G T (cr) 
is characterized by having i>_(er p ) = v-(cr) [ 106 1 . The fact that the smallest symplectic 
eigenvalue is the same for both partially transposed CMs entails that 

E T (<r) = G r (cr) = [max{0, - logi>_(cr)}] 2 . (107) 

We thus consistently retrieve for the Gaussian contangle (or, equivalently, the contangle, as 
they coincide in this specific case), the expression previously found for the mixed symmetric 
reductions of fully symmetric three-mode pure states, Eq. d99l ). 

To our aims, it is useful here to provide a compact, operative definition of the 
Gaussian contangle for 1 x JV bipartite Gaussian states, based on the evaluation of Gaussian 
entanglement measures in section 14.5.21 If tr^y is the CM of a (generally mixed) bipartite 
Gaussian state where subsystem Si comprises one mode only, then the Gaussian contangle 
G T can be computed as 

G T (cr lU ) = G T (cr°p = g\m\^ g[x] = arcsinh 2 [v^T] . (108) 

Here cr°j^* corresponds to a pure Gaussian state, and m^y = m(cr°y) = y^Det er° pt = 
^Det cr° pt , with er°^* being the reduced CM of subsystem Si (Sj) obtained by tracing 
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over the degrees of freedom of subsystem Sj (Si). The CM cr° p denotes the pure bipartite 
Gaussian state which minimizes m(er P y) among all pure-state CMs o - ?, . such that a P y < 

<Ji\j- If <Tj|j is a pure state, then cr°y = er^j, while for a mixed Gaussian state Eq. (1 1 08b 
is mathematically equivalent to constructing the Gaussian convex roof. For a separable 
state m(cr°y) = 1. Here we have implicitly used the fact that the smallest symplectic 
eigenvalue £>_ of the partial transpose of a pure 1 x N Gaussian state cr p ^ is equal to 

v_ = ^/DcToi — VDct Ci — 1> as follows by recalling that the 1 x JV entanglement is 
equivalent to a 1 x 1 entanglement by virtue of the phase-space Schmidt decomposition (see 
section [T.4.1b and by exploiting Eq. ( f80b with A = 2, /z = 1 and \x\ = pi = 1/ \/Det er;. 

The Gaussian contangle G T , like all members of the family of measures of entanglement 
(see section l4.5.2b is completely equivalent to the Gaussian entanglement of formation [ 106 1, 
which quantifies the cost of creating a given mixed Gaussian state out of an ensemble of pure, 
entangled Gaussian states. 

Gaussian tangle. — Analogously, for alxJV bipartition associated to a pure Gaussian 
state q p A \ B with Sa = S\ (a subsystem of a single mode) and Sb — <?2 • • ■ Sn, we define the 
following quantity 

tg(9 p aib )=M 2 (q p aib ). (109) 

Here, Af(g) is the negativity, Eq. d63l l, of the Gaussian state g. The functional tq, like the 
negativity TV, vanishes on separable states and does not increase under LOCC, i.e. , it is a 
proper measure of pure-state bipartite entanglement. It can be naturally extended to mixed 
Gaussian states pa\b via the convex roof construction 

tg(qa\b)= inf Vp l T G (pf), (110) 

{P„4 P) } i 

where the infimum is taken over all convex decompositions of gA\B m terms of pure Gaussian 
states g p : pa\b = Yli Vi(Fi - By virtue of the Gaussian convex roof construction, the Gaussian 
entanglement measure tq Eq. ( II 10b is an entanglement monotone under Gaussian LOCC (see 
section l4.5.21 i. Henceforth, given an arbitrary TV-mode Gaussian state fts 1 |,s 2 ...s JV , we refer to 
tq, Eq. ( 11 lOt as the Gaussian tangle |122|. Obviously, in terms of CMs, the analogous of 
the definition d 1 06b is valid for the Gaussian tangle as well, yielding it computable like the 
contangle in Eq. d 1 08b - Namely, exploiting Eq. d67l ), one finds 

1 _ 2 

TG((Ti\i)=TG((T°$)= l w[m\ j }, w[x] = - {\/ X - I + s/x - l) ,(111) 

where we refer to the discussion immediately after Eq. ( 1108b for the definition of the quantities 
involved in Eq. ( Ill lb . 

We will now proceed to review the structure of entanglement sharing in Gaussian 
states and the expressions of the monogamy constraints on its distribution. We remark that, 
being the (squared) negativity a monotonic and convex function of the (squared) logarithmic 
negativity, see Eq. ( 1102b . the validity of any monogamy constraint on distributed Gaussian 
entanglement using as a measure of entanglement the "Gaussian tangle", is implied by the 
proof of the corresponding monogamy inequality obtained using the "Gaussian contangle". 
For this reason, when possible, we will always employ as a preferred choice the primitive 
entanglement monotone, represented by the (Gaussian) contangle 111211 (which could be 
generally referred to as a 'logarithmic' tangle in quantum systems of arbitrary dimension). 
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7.3. Monogamy inequality for all Gaussian states. 

We are now in the position to recall a collection of recent results concerning the monogamy 
of distributed Gaussian entanglement in multimode Gaussian states. 

In the broadest setting we want to investigate whether a monogamy inequality like 
Ineq. d97| ) holds in the general case of Gaussian states with an arbitrary number TV of modes. 
Considering a Gaussian state distributed among N parties (each owning a single mode), the 
monogamy constraint on distributed entanglement can be written as 



where the global system is multipartitioned in subsystems (k = 1, . . ., TV), each owned by 
a respective party, and E is a proper measure of bipartite entanglement. The corresponding 
general monogamy inequality is known to hold for qubit systems 11 1 341 . 

The left-hand side of inequality ( II 12b quantifies the bipartite entanglement between a 
probe subsystem 5, and the remaining subsystems taken as a whole. The right-hand side 
quantifies the total bipartite entanglement between 5, and each one of the other subsystems 
Sjjti in the respective reduced states. The non negative difference between these two 
entanglements, minimized over all choices of the probe subsystem, is referred to as the 
residual multipartite entanglement. It quantifies the purely quantum correlations that are not 
encoded in pairwise form, so it includes all manifestations of genuine if -partite entanglement, 
involving K subsystems (modes) at a time, with 2 < K < N. The study of entanglement 
sharing and monogamy constraints thus offers a natural framework to interpret and quantify 
entanglement in multipartite quantum systems [ 136|. 

With these premises, we have proven that the (Gaussian) contangle (and the Gaussian 
tangle, as an implication) is monogamous in fully symmetric Gaussian states of N modes 
11211 . In general, we have proven the Gaussian tangle to satisfy inequality (11 121 i in all 
TV-mode Gaussian states 01221 . A full analytical proof of the monogamy inequality for 
the contangle in all Gaussian states beyond the symmetry, is currently lacking; however, 
numerical evidence obtained for randomly generated non symmetric 4-mode Gaussian states 
strongly supports the conjecture that the monogamy inequality be true for all multimode 
Gaussian states, using also the (Gaussian) contangle as a measure of bipartite entanglement 
11211 . Remarkably, for all (generally nonsymmetric) three-mode Gaussian states the 
(Gaussian) contangle has been proven to be monogamous, leading in particular to a proper 
measure of tripartite entanglement in terms of residual contangle: the analysis of distributed 
entanglement in the special instance of three-mode Gaussian states, with all the resulting 
implications, is postponed to the next section. 

Let us restate again the main result of this section: The Gaussian tangle tq, an 
entanglement monotone under Gaussian LOCC, is monogamous for all, pure and mixed, N- 
mode Gaussian states distributed among N parties, each owning a single mode A122I . 

7.4. Discussion. 

The monogamy constraints on entanglement sharing are essential for the security of CV 
quantum cryptographic schemes 11381 11391 . because they limit the information that might 
be extracted from the secret key by a malicious eavesdropper. Monogamy is useful as well 
in investigating the range of correlations in Gaussian matrix-product states of harmonic rings 
Ml 401 l30l . and in understanding the entanglement frustration occurring in ground states of 



JY 




(112) 
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many-body harmonic lattice systems 11411 . which may be now extended to arbitrary states 
beyond symmetry constraints. 

At a fundamental level, having established the monogamy property for all Gaussian states 
paves the way to a proper quantification of genuine multipartite entanglement in CV systems 
in terms of the residual distributed entanglement. In this respect, the intriguing question arises 
whether a stronger monogamy constraint exists on the distribution of entanglement in many- 
body systems, which imposes a physical trade-off on the sharing of both bipartite and genuine 
multipartite quantum correlations. 

It would be important to understand whether the inequality dl 121 ) holds as well for 
discrete-variable qudits (2 < d < oo), interpolating between qubits and CV systems 111361 . 
If this were the case, the (convex-roof extended) squared negativity, which coincides with 
the tangle for arbitrary states of qubits and with the Gaussian tangle for Gaussian states of 
CV systems, would qualify as a universal bona fide, dimension-independent quantifier of 
entanglement sharing in all multipartite quantum systems. At present, this must be considered 
as a completely open problem. 

8. Multipartite Gaussian entanglement: quantification and promiscuous sharing 
structure. 

This section is mainly devoted to the characterization of entanglement in the simplest 
multipartite CV setting, namely that of three-mode Gaussian states. 

To begin with, let us set the notation and review the known results about three-mode 
Gaussian states of CV systems. We will refer to the three modes under exam as mode 1, 2 
and 3. The 2x2 submatrices that form the CM er = er 12 3 of a three-mode Gaussian state are 
defined according to Eq. ( fT9l l. whereas the 4x4 CMs of the reduced two-mode Gaussian states 
of modes i and j will be denoted by er y . Likewise, the local (two-mode) seralian invariants 
Ay, Eq. d32l . will be specified by the labels i and j of the modes they refer to, while, to 
avoid any confusion, the three-mode (global) seralian symplectic invariant will be denoted by 
A = A123. Let us recall the uncertainty relation Eq. d73l for two-mode Gaussian states, 

A y - Deter y < 1 . (113) 
8.1. Separability classes for three-mode Gaussian states. 

As it is clear from the discussion of section 14.4.11 a complete qualitative characterization 
of the entanglement of three-mode Gaussian state is possible because the PPT criterion is 
necessary and sufficient for their separability under any, partial or global (i.e. 1 x 1 or 1 x 2), 
bipartition of the modes. This has lead to an exhaustive classification of three-mode Gaussian 
states in five distinct separability classes II 1421 . These classes take into account the fact that 
the modes 1, 2 and 3 allow for three distinct global bipartitions: 

• Class 1: states not separable under all the three possible bipartitions i x (Jk) of the modes 
(fully inseparable states, possessing genuine multipartite entanglement). 

• Class 2: states separable under only one of the three possible bipartitions (one-mode 
biseparable states). 

• Class 3: states separable under only two of the three possible bipartitions (two-mode 
biseparable states). 

• Class 4: states separable under all the three possible bipartitions, but impossible to write 
as a convex sum of tripartite products of pure one-mode states (three-mode biseparable 
states). 
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• Class 5: states that are separable under all the three possible bipartitions, and can be 
written as a convex sum of tripartite products of pure one-mode states (fully separable 
states). 

Notice that Classes 4 and 5 cannot be distinguished by partial transposition of any of the 
three modes (which is positive for both classes). States in Class 4 stand therefore as 
nontrivial examples of tripartite entangled states of CV systems with positive partial transpose 
[ 1 42 1 . It is well known that entangled states with positive partial transpose possess bound 
entanglement, that is, entanglement that cannot be distilled by means of LOCC. 



8.2. Residual contangle as genuine tripartite entanglement monotone. 

We have proven in Ref. 111211 . that all three-mode Gaussian states satisfy the CKW monogamy 
inequality d97l i. using the (Gaussian) contangle Eq. ( 1108b to quantify bipartite entanglement. 

The sharing constraint leads naturally to the definition of the residual contangle as a 
quantifier of genuine tripartite entanglement in three-mode Gaussian states, much in the same 
way as in systems of three qubits fl 1 3 3H . However, at variance with the three-qubit case (where 
the residual tangle of pure states is invariant under qubit permutations), here the residual 
contangle is partition-dependent according to the choice of the probe mode, with the obvious 
exception of the fully symmetric states. A bona fide quantification of tripartite entanglement 
is then provided by the minimum residual contangle 111 211 



E$ j \ k = min 



(114) 



where the symbol (i,j,k) denotes all the permutations of the three mode indexes. This 
definition ensures that E l }^ k is invariant under all permutations of the modes and is thus a 
genuine three-way property of any three-mode Gaussian state. We can adopt an analogous 
definition for the minimum residual Gaussian contangle G r T es , sometimes referred to as 
arravogliament II121II1311I1431 : 



G res _ G i\j\k = mm Qi\(jk) _ Qi\j _ QAk 



(115) 



One can verify that 

^Qi\{jk) _ Qijfcj _ (Qj\{ik) _ G j|fc) > (U6) 

if and only if > aj, and therefore the absolute minimum in Eq. ( II 14b is attained by the 
decomposition realized with respect to the reference mode I of smallest local mixedness ai, 
i.e. for the single-mode reduced state with CM of smallest determinant (corresponding to the 
largest local purity /!;)■ 

A crucial requirement for the residual (Gaussian) contangle, Eq. ( II 15b , to be a proper 
measure of tripartite entanglement is that it be nonincreasing under (Gaussian) LOCC. The 
monotonicity of the residual tangle was proven for three-qubit pure states in Ref. 111441 . 
In the CV setting we have proven that for pure three-mode Gaussian states G r T es is an 
entanglement monotone under tripartite Gaussian LOCC, and that it is nonincreasing even 
under probabilistic operations, which is a stronger property than being only monotone on 
average 11211 . Therefore, the residual Gaussian contangle G r T es is a proper and computable 
measure of genuine multipartite (specifically, tripartite) entanglement in three-mode Gaussian 
states. 
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8.3. Standard form and tripartite entanglement of pure three-mode Gaussian states. 

Here we apply the above defined measure of tripartite entanglement to the relevant instance 
of pure three-mode Gaussian states. We begin by recalling their structural properties. 

8.3.1. Symplectic properties of three-mode covariance matrices. The CM <x of a pure three- 
mode Gaussian state is characterized by 

Det<r = l, A = 3. (117) 

The purity constraint requires the local entropic measures of any 1 x 2-mode bipartitions to 
be equal: 

Deter^ = Deter fe , (118) 

with i, j and k different from each other. This general, well known property of the bipartitions 
of pure states may be easily proven resorting to the Schmidt decomposition (see section [T.4.1| i. 

A first consequence of Eqs. dl 171 ) and dl 181 ) is rather remarkable. Combining such 
equations one easily obtains 

(Aia - Det er 12 ) + (A 13 - Det er 13 ) + (A 23 - Det <x 23 ) = 3 , 

which, together with Inequality dl 131 ), implies 

A y = Det (Tij + 1 , Vt, j : i ^ ] . (119) 

The last equation shows that any reduced two-mode state of a pure three-mode Gaussian 
state saturates the partial uncertainty relation Eq. dl 131 ). The states endowed with such a 
partial minimal uncertainty (namely, with their smallest symplectic eigenvalue equal to 1) are 
states of minimal negativity for given global and local purities, alias GLEMS (Gaussian least 
entangled mixed states) l|104|[3"71 , introduced in section [53] 

This is relevant considering that the standard form CM of Gaussian states is completely 
determined by their global and local invariants. Therefore, because of Eq. dl 18t , the 
entanglement between any pair of modes embedded in a three-mode pure Gaussian state is 
fully determined by the local invariants Dot er/, for I = 1,2,3, whatever proper measure we 
choose to quantify it. Furthermore, the entanglement of a <j; L \<jjk bipartition of a pure three- 
mode state is determined by the entropy of one of the reduced states that is, once again, by 
the quantity Det cr,;. Thus, the three local symplectic invariants Det cr\, Det er 2 and Det er 3 
fully determine the entanglement of any bipartition of a pure three-mode Gaussian state. We 
will show that they suffice to determine as well the genuine tripartite entanglement encoded 
in the state II 1 3 1 H . 

For ease of notation, in the following we will denote by a/ the local single-mode 
symplectic eigenvalues associated to mode I with CM erf. 

ai = VDct en . (120) 

Eq. (|46| l shows that the quantities a; are simply related to the purities of the reduced single- 
mode states, the local purities fii, by the relation 

1*1 = -- (121) 

ai 

Since the set {a{\ fully determines the entanglement of any of the 1x2 and lxl bipartitions 
of the state, it is important to determine the range of the allowed values for such quantities. 
This will provide a complete quantitative characterization of the entanglement of three-mode 
pure Gaussian states. 
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Such an analysis has been performed in Ref . [131]. Defining the parameters 

a[ = a t - 1 , (122) 
one has that the following permutation-invariant triangular inequality 

K-af 3 \ < af k < a^ + af J . (123) 

holds as the only condition, together with the positivity of each aj, which fully characterizes 
the local symplectic eigenvalues of the CM of three-mode pure Gaussian states. All standard 
forms of pure three-mode Gaussian states and in particular, remarkably, all the possible values 
of the negativities (section 14.5. Q and/or of the Gaussian entanglement measures (section 
14.5.2b between any pair of subsystems, can be determined by letting a[, a' 2 and a' 3 vary in 
their range of allowed values. Let us remark that Eq. ( I123l l qualifies itself as an entropic 
inequality, as the quantities {a'j} are closely related to the purities and to the von Neumann 
entropies of the single-mode reduced states. In particular the von Neumann entropies Svj of 
the reduced states are given by Syj = f(a'j + 1) = f{o-j), where the increasing convex 
entropic function f{x) has been defined in Eq. d481 l. Now, Inequality (11231 > is strikingly 
analogous to the well known triangle (Araki-Lieb) and subadditivity inequalities for the von 
Neumann entropy (holding for general systems [ 145 1 36)), which in our case read 

|/(0i) - /( aj )| < f(a k ) < f(cn) + f{aj) . (124) 

However, the condition imposed by Eq. ( I123l l is strictly stronger than the generally holding 
inequalities ( I124l i for the von Neumann entropy applied to pure quantum states. 

We recall that the form of the CM of any Gaussian state can be simplified through local 
(unitary) symplectic operations, that therefore do not affect the entanglement or mixedness 
properties of the state, belonging to Spf^ |f30l, 146, 106|. Such reductions of the CMs are 
called "standard forms", as introduced in section [3~4l For the sake of clarity, let us write the 
explicit standard form CM of a generic pure three-mode Gaussian state II 1 3 11 , 
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with 



(Oj - a/) 2 - (a k - l) 2 



(di - dj) 2 - (a k + l) 2 



(126) 



± 



) 2 - (ofc - l) 2 



(oj + aj) 2 - (a k + l) 2 



Let us stress that, although useful in actual calculations, the use of CMs in standard form 
does not entail any loss of generality, because all the results concerning the characterization 
of bipartite and tripartite entanglement do not depend on the choice of the specific form of the 
CMs, but only on invariant quantities, such as the global and local symplectic invariants. 
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8.3.2. Residual contangle of fully inseparable three-mode Gaussian states. We can now 
analyze separability and entanglement in more detail. From our study it turns out that, 
regarding the classification of section [84] I 142], pure three-mode Gaussian states may belong 
either to Class 5, in which case they reduce to the global three-mode vacuum, or to Class 2, 
reducing to the uncorrelated product of a single-mode vacuum and of a two-mode squeezed 
state, or to Class 1 (fully inseparable state). No two-mode or three-mode biseparable pure 
three-mode Gaussian states are allowed. 

Let us now describe the complete procedure to determine the genuine tripartite 
entanglement in a pure three-mode Gaussian state with a completely general (not necessarily 
in standard form) CM <t p belonging to Class 1, as presented in Ref. 11311 . 

(i) Determine the local purities The state is globally pure (Det cr p = 1). The only quantities 

needed for the computation of the tripartite entanglement are therefore the three local 
mixednesses a;, defined by Eq. dl20t , of the single-mode reduced states <r ; , / = 1, 2, 3 
(see Eq. (fT9l). Notice that the global CM cr p needs not to be in the standard form 
of Eq. (11251 ). as the single-mode determinants are local symplectic invariants. From 
an experimental point of view, the parameters a; can be extracted from the CM using 
the homodyne tomographic reconstruction of the state 111 171 : or they can be directly 
measured with the aid of single photon detectors 111 II II 1211 . 

(ii) Find the minimum From Eq. ( 11 161 ), the minimum in the definition dl 151 ) of the residual 

Gaussian contangle G r T es is attained in the partition where the bipartite entanglements 
are decomposed choosing as probe mode I the one in the single-mode reduced state of 
smallest local mixedness a; = a m i n . 

(iii) Check range and compute Given the mode with smallest local mixedness a m i n (say, for 
instance, mode 1) and the parameters s and d defined by 

0,2 + a.3 



h° Q>min 

then 



a-2 - a 3 



1 then mode 1 is uncorrelated from the others: G 1 ! 



(127) 
(128) 

0. If, instead, a m i n > 1 



with Q = Gt 2 



G r T es (crP) = arcsinh' 
Gt' 3 defined by 



a min 1 



Q(a, s, d) — arcsinh 2 y m 2 (a, s, d) — 1 
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s/ m 2 ( 



a, s, 



-d) - 1 



where m = to_ if D < 0, and m = m + otherwise (one has 
Here: 

|fc-l 

m - (s - df - 1 ' 



m_ for D = 0). 



2a 2 (l + 2s 2 + 2d 2 ) - (4s 2 - l)(4d 2 - 1) - a 4 - 
4(s - d) 



D = 2{s - d) - ^2 [fc 2 + 2k+ + \k-\{k 2 _ + 8fc+) 1 /2] /k+ 
k ± =a 2 ±(s + d) 2 , 




Figure 7. Three-dimensional plot of the residual Gaussian contangle G T es (<r p ) in pure three- 
mode Gaussian states rr p , determined by the three local mixednesses a;, I = 1, 2, 3. One of 
the local mixednesses is kept fixed (ai = 2). The remaining ones vary constrained by the 
triangle inequality )123t . The explicit expression of G r es is given by Eq. 4 1 29i . See text for 
further details. 



5 =(a- 2d-l)(a- 2d + l)(a + 2d - l)(a + 2d+l) 
(a - 2s - l)(a - 2s + l)(a + 2s - l)(a + 2s + 1) . 

Note (we omitted the explicit dependence for brevity) that each quantity in Eq. d 1 3 1 b 
is a function of (a, s, d). Therefore, to evaluate the second term in Eq. ( 1130b each d in 
Eq. ( 1131b must be replaced by — d. Note also that if d < — {af nm — l) /As then G r ' 2 = 0. 
Instead, if d > {a^ nin — l)/4s then G r ' 3 = 0. Otherwise, all terms in Eq. ( II 15b are 
nonvanishing. 

The residual Gaussian contangle Eq. ( II 151 ) in generic pure three-mode Gaussian states is 
plotted in Fig.|7]as a function of 02 and 03, at constant a\ — 2. For fixed a\, it is interesting 
to notice that G r es is maximal for a% = 0,3, i.e. for bisymmetric states. Notice also how the 
residual Gaussian contangle of these bisymmetric pure states has a cusp for a\ = 0,2 — a.3. 
In fact, from Eq. dl 16b , for 02 = 03 < ai the minimum in Eq. dl 151 ) is attained decomposing 
with respect to one of the two modes 2 or 3 (the result is the same by symmetry), while for 
0-2 = 0-3 > a,\ mode 1 becomes the probe mode. 

8.3.3. Residual contangle and distillability of mixed states. For generic mixed three-mode 
Gaussian states, a quite cumbersome analytical expression for the 1|2 and 1 1 3 Gaussian 
contangles may be written, which explicitly solves the minimization over the angle 9 in 
Eq. d90b . On the other hand, the optimization appearing in the computation of the 1 1 (23) 
bipartite Gaussian contangle (see Eq. ( 1108b ) has to be solved only numerically. However, 
exploiting techniques like the unitary localization of entanglement described in section [6] 
and results like that of Eq. ( 1107b . closed expressions for the residual Gaussian contangle 
can be found as well in relevant classes of mixed three-mode Gaussian states endowed with 
some symmetry constraints. Interesting examples of these states and the investigation of their 
physical properties are discussed in Refs. Ill31l fl43l. 
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As an additional remark, let us recall that, although the entanglement of Gaussian states is 
always distillable with respect to 1 x JV bipartitions [ 80 1 (see section |4~.4.1| ), they can exhibit 
bound entanglement in 1 x 1 x 1 tripartitions 111421 . In this case, the residual Gaussian 
contangle cannot detect tripartite PPT entangled states. For example, the residual Gaussian 
contangle in three-mode biseparable Gaussian states (Class 4 of Ref. B142I ) is always zero, 
because those bound entangled states are separable with respect to all 1 x 2 bipartitions of the 
modes. In this sense we can correctly regard the residual Gaussian contangle as an estimator 
of distillable tripartite entanglement, being strictly nonzero only on fully inseparable three- 
mode Gaussian states (Class 1 in the classification of section [8~TT >. 



8.4. Sharing structure of multipartite entanglement: promiscuous Gaussian states. 

We are now in the position to review the sharing structure of CV entanglement in three- 
mode Gaussian states by taking the residual Gaussian contangle as a measure of tripartite 
entanglement, in analogy with the study done for three qubits [ 144 1 using the residual tangle 

cm 

The first task we face is that of identifying the three-mode analogues of the two 
inequivalent classes of fully inseparable three-qubit states, the GHZ state II 1471 

I^ghz) = ^= (|000) + , (131) 

and the W state 01441 

\^ w ) = (|001) + |010) + |100)) . (132) 

These states are both pure and fully symmetric, i.e. invariant under the exchange of any two 
qubits. On the one hand, the GHZ state possesses maximal tripartite entanglement, quantified 
by the residual tangle 111331 [144], with zero couplewise entanglement in any reduced state 
of two qubits reductions. Therefore its entanglement is very fragile against the loss of 
one or more subsystems. On the other hand, the W state contains the maximal two-party 
entanglement in any reduced state of two qubits [ 144] and is thus maximally robust against 
decoherence, while its tripartite residual tangle vanished 



8.4.1. CV finite-squeezing GHZ/SN states. To define the CV counterparts of the three-qubit 
states |^ghz) an d |'0w)> one must start from the fully symmetric (generally mixed) three- 
mode CM <j s of the form cr a 3, Eq. ( Bil l. Surprisingly enough, in symmetric three-mode 
Gaussian states, if one aims at maximizing, at given single-mode mixedness a = VDct a, 
either the bipartite entanglement Gffi in any two-mode reduced state (i.e. aiming at the CV W- 
like state), or the genuine tripartite entanglement G r T es {i.e. aiming at the CV GHZ-like state), 
one finds the same, unique family of states. They are exactly the pure, fully symmetric three- 
mode Gaussian states (three-mode squeezed states) with CM cr p s of the form cr a 3, Eq. ( BTT i. 
with a = al 2 , e = diag{e + , e~} and 



, a 2 - 1 ± J(a 2 - 1) (9a 2 - 1) 

e ± = ^ '- , (133) 

4a 

where we have used Eq. d94i i ensuring the global purity of the state. In general, we have 
studied the entanglement scaling in fully symmetric (pure) Gaussian states by means of 
the unitary localization in section 16.21 It is in order to mention that these states were 

+ The full inseparability of the W state can be however detected by the 'Schmidt measure' 1 148] or by introducing 
also for qubits a monogamy inequality in terms of negativity 11491 . 
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previously known in the literature as CV "GHZ-type" states [128 90], as in the limit of 
infinite squeezing (a — > oo), they approach the proper (unnormalizable) continuous-variable 
GHZ state J dx\x,x,x), a simultaneous eigenstate of total momentum pi + p 2 + pz and of 
all relative positions q\ — q~j = 1, 2, 3), with zero eigenvalues [ 1 50 1 . 

For any finite squeezing (equivalently, any finite local mixedness a), however, the above 
entanglement sharing study leads ourselves to re-baptize these states as "CV GHZ/W states" 
EI] EI] El, and denote their CM by aT z/w ■ 

The residual Gaussian contangle of GHZ/W states of finite squeezing takes the simple 
form ED 



G; es (<T? HZ/B ') = arcsmh 2 y/tf^l 



(134) 



3a 2 - 1 - V9a 4 - 10a 2 + 1 



It is straightforward to see that G£ es (<Tg w ) is nonvanishing as soon as a > 1. 
Therefore, the GHZ/W states belong to the class of fully inseparable three-mode states 
[142] [128] ED HQ) (Class 1, see section ED. We finally recall that in a GHZ/W state the 
residual Gaussian contangle G r T es Eq. d 1 15b coincides with the true residual contangle E^' 2 ' 3 
Eq. (II 141 . This property clearly holds because the Gaussian pure-state decomposition is the 
optimal one in every bipartition, due to the fact that the global three-mode state is pure and 
the reduced two-mode states are symmetric (see section l5".2.2b . 

The peculiar nature of entanglement sharing in CV GHZ/W states is further confirmed 
by the following observation. If one requires maximization of the 1x2 bipartite Gaussian 

i ( jk) 

contangle G r under the constraint of separability of all the reduced two-mode states, one 
finds a class of symmetric mixed states (T states 112111131111431 ') whose residual contangle 
is strictly smaller than the one of the GHZ/W states Eq. ( I134l i for any fixed value of the local 
mixedness a, that is, for any fixed value of the only parameter (operationally related to the 
squeezing of each single mode) that completely determines the CMs of both families of states 
up to local unitary operations. 



8.4.2. Promiscuous continuous-variable entanglement sharing. The above results lead to 
the conclusion that in symmetric three-mode Gaussian states, when there is no bipartite 
entanglement in the two-mode reduced states (like in T states) the genuine tripartite 
entanglement is not enhanced, but frustrated. More than that, if there are maximal quantum 
correlations in a three-party relation, like in GHZ/W states, then the two-mode reduced states 
of any pair of modes are maximally entangled mixed states. 

These findings establish the promiscuous nature of CV entanglement sharing in 
symmetric Gaussian states 11211 . Being associated with degrees of freedom with continuous 
spectra, states of CV systems need not saturate the CKW inequality to achieve maximum 
couplewise correlations (as it was instead the case for W states of qubits, Eq. ( 1132b ). In fact, 
without violating the monogamy constraint Ineq. (I97t . pure symmetric three-mode Gaussian 
states are maximally three-way entangled and, at the same time, possess the maximum 
possible entanglement between any pair of modes in the corresponding two-mode reduced 
states. The notion of "promiscuity" basically means that bipartite and genuine multipartite 
(in this case tripartite) entanglement are increasing functions of each other, while typically in 
low-dimensional systems like qubits only the opposite behavior is compatible with monogamy 
II 1361 . The promiscuity of entanglement in three-mode GHZ/W states is, however, partial. 
Namely they exhibit, with increasing squeezing, unlimited tripartite entanglement (diverging 
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in the limit a — > oo) and nonzero, accordingly increasing bipartite entanglement between 
any two modes, which nevertheless stays finite even for infinite squeezing. Precisely, from 
Eq. d!34t , it saturates to the value 

rii\ j I ghz/w \ log 3 

G T M (er s , a — > oo) = — - — f» 0.3 . (135) 

In Gaussian states of CV systems with more than three modes, entanglement can indeed be 
distributed in an infinitely promiscuous way 11321 . as we will briefly discuss in the following. 

More remarks are in order concerning the tripartite case. The structure of entanglement 
in GHZ/W states is such that, while being maximally three-party entangled, they are also 
maximally robust against the loss of one of the modes. This preselects GHZ/W states 
also as optimal candidates for carrying quantum information through a lossy channel, being 
intrinsically less sensitive to decoherence effects. They have been in fact proven to be 
maximally robust against environmental noise among all three-mode Gaussian states 01311 . 

It is natural to question whether all three-mode Gaussian states are expected to exhibit a 
promiscuous entanglement sharing. Such a question is addressed in [ 143 1, by investigating the 
persistency of promiscuity against the lack of each of the two defining properties of GHZ/W 
states: full symmetry, and global purity. One specifically finds that entanglement promiscuity 
survives under a quite strong amount of mixedness (up to an impurity of 1 — \i « 0.8, 
but is in general lost if the complete permutation-invariance is relaxed. Pure three-mode 
Gaussian states which are only bisymmetric and not fully symmetric (known as basset hound 
states) offer indeed a traditional, not promiscuous entanglement sharing, with bipartite and 
tripartite entanglement being competitors 01431 . Therefore, in the tripartite Gaussian setting, 
'promiscuity' is a peculiar consequence not of the global purity, but of the complete symmetry 
under modes-exchange. Beside frustrating the maximal entanglement between pairs of modes 
II 1411 . symmetry also constrains the multipartite sharing of quantum correlations. In basset 
hound states, the separability of the reduced state of modes 2 and 3 prevents the three modes 
from having a strong genuine tripartite entanglement among them all, despite the heavy 
quantum correlations shared by the two couples of modes 1\2 and 1|3. 

8.4.3. Unlimited promiscuity of entanglement in four-mode Gaussian states. The above 
argument on the origin of promiscuity does not hold in the case of Gaussian states with four 
and more modes, where relaxing the symmetry constraints may allow for an enhancement 
of the distributed entanglement promiscuity to an unlimited extent. In Ref. [132|, we have 
introduced a class of pure four-mode Gaussian states which are not fully symmetric, but 
invariant under the double exchange of modes 1 <-> 4 and 2^3. They are defined as 
follows. 

One starts with an uncorrelated state of four modes, each one initially in the vacuum 
of the respective Fock space, whose corresponding CM is the identity. One applies a two- 
mode squeezing transformation 6*2,3(5), Eq. (1231 . with squeezing s to modes 2 and 3, then 
two further two-mode squeezing transformations Si^(ct) and 53.4(a), with squeezing a, to 
the pairs of modes {1,2} and {3,4}. The two last transformations serve the purpose of 
redistributing the original bipartite entanglement, created between modes 2 and 3 by the first 
two-mode squeezing operations, among all the four modes. For any value of the parameters s 
and a, the output is a pure four-mode Gaussian state with CM cr, 

cr = S 3 A{a)S ly 2(a)S 2 ,3(s)S2 i3 (s)Sj i2 (a)Sj A (a) . (136) 

Explicitly, cr is of the form Eq. ( fT~9b where 

<j\ = er 4 = [cosh 2 (a) + cosh(2s) sinh 2 (a)]l2 , 
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<r 2 = o"3 = [cosh(2s) cosh 2 (a) + sinh 2 (a)]l2 , 

£1,2 = £3,4 = [cosh 2 (s) sinh(2a)]Z 2 , 

£1,3 = £2,4 = [cosh(a) sinh(a) smh(2s)]l 2 , 

£l t 4 = [sinh 2 (a) sinh(2s)]Z 2 , 

e 2j 3 = [cosh 2 (a) sinh(2s)]Z 2 , 

with^^o-i)' 

By a proper investigation of the properties of distributed entanglement in the parametric 
state Eq. ( 11361 l. it can be shown that the entanglement between modes 1 and 2 is an 
unboundedly increasing function of the squeezing a; the same holds for the pair of modes 
3 and 4. In addition, the block of modes (1,2) is arbitrarily entangled with the block of modes 
(3,4) as a function of s. On the other hand, one can demonstrate that in the same state an 
unlimited genuine four-partite entanglement is present, increasing arbitrarily with increasing 
a, and thus coexisting with (and being mutually enhanced by) the bipartite entanglement in 
two pairs of modes 11321 . 

Such a simple example demonstrates that, when the quantum correlations arise among 
degrees of freedom spanning an infinite-dimensional space of states (characterized by 
unbounded mean energy), an accordingly infinite freedom is allowed for different forms of 
bipartite and multipartite entanglement. This phenomenon happens with no violation of the 
fundamental monogamy constraint that retains its general validity in quantum mechanics. In 
the CV instance the only effect of monogamy is to bound the divergence rates of the individual 
entanglement contributions as the squeezing parameters are increased. Within the restricted 
Hilbert space of four or more qubits, instead, an analogous entanglement structure between 
the single qubits is strictly forbidden. 

This result opens interesting perspectives for the understanding and characterization of 
entanglement in multiparticle systems. Gaussian states with finite squeezing (finite mean 
energy) are somehow analogous to discrete systems with an effective dimension related to the 
amount of squeezing [ 1 1. As the promiscuous entanglement sharing arises in Gaussian states 
by asymptotically increasing the squeezing to infinity, it is natural to expect that dimension- 
dependent families of states will exhibit an entanglement structure that becomes gradually 
more promiscuous with increasing dimension of the Hilbert space towards the CV limit. A 
proper investigation on systems of qudits (2 < d < 00) is therefore a necessary step in order 
to develop a complete picture of entanglement sharing in many-body systems [136|. This 
program has been initiated by establishing a sharp discrepancy between the two extrema in the 
ladder of Hilbert space dimensions: in the case of C V systems in the limit of infinite squeezing 
(infinite mean energy) entanglement has been proven infinitely more shareable than that of 
individual qubits Ml 321 . This fact could prelude to implementations of quantum information 
protocols with CV systems that cannot be achieved and not even devised with qubit resources. 

9. Conclusions and outlook. 

9.1. Entanglement in non Gaussian states. 

The infinite-dimensional quantum world is obviously not confined to Gaussian states. In fact, 
some recent results demonstrate that the current trends in the theoretical understanding and 
experimental control of CV entanglement are strongly pushing towards the boundaries of the 
territory of Gaussian states and Gaussian operations. 

The entanglement of Gaussian states cannot be increased (distilled) by Gaussian 
operations (83] [84] 1851 . Similarly, for universal one-way quantum computation using 
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Gaussian cluster states, a single-mode non Gaussian measurement is required 111521 . 
Moreover, a fundamental motivation for investigating entanglement in non Gaussian states 
stems from the property of extremality of Gaussian states: it has been recently proved that 
they are the least entangled states among all states of CV systems with given second moments 
1 1 5 3 1 . Experimentally, it has been recently demonstrated 015411 that a two-mode squeezed 
Gaussian state can be "degaussified" by coherent subtraction of a single photon, resulting in 
a mixed non Gaussian state whose nonlocal properties and entanglement degree are enhanced 
(enabling a better efficiency for teleporting coherent states I1155ID . Theoretically, even the 
characterization of bipartite entanglement (let alone multipartite) in non Gaussian states stands 
as a formidable task. 

One immediate observation is that any two-mode state with a CM corresponding to an 
entangled Gaussian state is itself entangled H 15 IB . Therefore, most of the results reviewed in 
this paper may serve to detect entanglement in a broader class of states of infinite-dimensional 
Hilbert spaces. They are, however, all sufficient conditions on entanglement based only on 
the second moments of the canonical operators. As such, for arbitrary non Gaussian states, 
they are in general very inefficient (meaning that most entangled non Gaussian states fail 
to be detected by these criteria). The description of non Gaussian states requires indeed to 
consider high order statistical moments. It could then be expected that inseparability criteria 
for these states should involve high order correlations. Recently, some separability criteria 
based on hierarchies of conditions involving higher moments of the canonical operators have 
been introduced to provide a sharper detection of inseparability in generic non Gaussian states. 

A first step in this direction has been taken by Agarwal and Biswas, who have applied 
the method of partial transposition to the uncertainty relations for the Schwinger realizations 
of the SU (2) and SU (1, 1) algebras 111561 . This approach can be successfully used to detect 
entanglement in the two-mode entangled non Gaussian state described by the wave function 
ip(x 1 ,x 2 ) = (2/7r) 1/2 (7ia;i + 7 2 x 2 )ea;p-(a:f + a;|)/2, with |7i| 2 + I72I 2 = 1, and in the 
SU(2) minimum-uncertainty states 11561 11571 . The demonstration of criteria consisting in 
hierarchies of arbitrary order moments has been achieved preliminarily by Hillery and Zubairy 
11581 , and definitively by Shchukin and Vogel 01591 11601 1161| . In particular, Shchukin and 
Vogel have introduced an elegant and unifying approach, based on the PPT requirement, that 
allows to derive, in the form of an infinite series of inequalities, a necessary and sufficient 
condition for the negativity of the partial transposition q of a bipartite quantum state q. The 
Shchukin- Vogel (SV) criterion includes as special cases all the above mentioned conditions 
(including the ones on second moments ll79l[T3l qualifying entanglement in Gaussian states), 
thus demonstrating the important role of PPT in building a strong criterion for the detection 
of entanglement. 

Here we briefly review the main features of the SV criterion first introduced by 
Shchukin and Vogel II 1591 . and later analyzed in mathematical detail by Miranowicz and 
Piani 1 160, 161] (see also 1 162|). Consider two bosonic modes A\ and A2 with the associated 
annihilation and creation operators, respectively ai,a\ and 0,2,0,2^- Shchukin and Vogel 
showed that every Hermitian operator X that acts on the Hilbert space of the two modes 
is nonnegative if and only if for any operator / whose normally-ordered form exists, i.e., for 
any operator / that can be written in the form 




(137) 



it holds that 



Tr/t/i" > 0. 



(138) 
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Applying the above result to the partially transposed matrix p of a quantum state density 
matrix p, one finds that p is nonpositive if and only if 

Tr[e/t/]<0. (139) 

Eq. ( 11391 ) is thus a necessary and sufficient condition for the negativity of the partially 
transposed density matrix (NPT). As such, the NPT criterion is in principle able to detect the 
bipartite entanglement of all CV two-mode states, pure or mixed, Gaussian or non Gaussian, 
with the exception of the PPT bound entangled states. 

Some specific applications to selected classes of non Gaussian states have been recently 
discussed in the two-mode setting [ 163 1, and preliminary extensions of the SV NPT criterion 
to multimode and multipartite cases have been introduced and applied to distinguish between 
different classes of separability B164II . To this aim, entanglement witnesses are useful as well 
||89l . Another, inequivalent condition based on matrices of moments (realignment criterion) 
has also been recently introduced for the qualification of bipartite entanglement in non 
Gaussian states (and in general mixed states of quantum systems in arbitrary dimension) II162I . 
We should mention a further interesting approach to non Gaussian entanglement analyzed by 
McHugh et al. 111651 , who showed that the entanglement of multiphoton squeezed states is 
completely characterized by observing that with respect to a new set of modes, those non 
Gaussian states actually assume a Gaussian character. 

The efficiency of some of the above-mentioned inseparability criteria in detecting the 
entanglement of non Gaussian, squeezed number states of two-mode radiation fields has been 
recently evaluated |166|. Detailed studies of the entanglement properties of non Gaussian 
states associated to SU(1, 1) active and SU(2) passive optical transformations, and the 
efficiency of the SV NPT criterion to reveal them are currently under way, together with 
the analysis of their dynamical properties in the presence of noise and decoherence 11671 . 

9.2. Applications, open problems, and current perspectives. 

The centrality of Gaussian states in CV quantum information is motivated not only by their 
peculiar structural properties which make their description amenable of an analytical analysis, 
but also by the ability to produce, manipulate and detect such states with remarkable accuracy 
in realistic, experimental settings. 

The scope of this review has been almost entirely theoretical. For reasons of space and 
time, we cannot discuss in sufficient detail all the proposals and experimental demonstrations 
concerning on one hand the state engineering of two-, three- and in general A-mode Gaussian 
states, and on the other hand the use of such states as resources for the realization of quantum 
information protocols. Excellent review papers are already available for what concerns both 
the optical state engineering of Gaussian and non Gaussian quantum states of CV systems 
1 168|, and the implementations of quantum information and communication with continuous 
variables 11511 [T). We just mention that, concerning the state engineering, an efficient 
scheme to produce generic pure A-mode Gaussian states in a standard form not encoding 
direct correlations between position and momentum operators (and so encompassing all the 
instances of multimode Gaussian states introduced in the previous sections) has been recently 
proposed 1 30 ] ; it enables to interpret entanglement in this subclass of Gaussian states entirely 
in terms of the two-point correlations between any pair of modes. 

From a practical point of view, Gaussian resources have been widely used to implement 
paradigmatic protocols of CV quantum information, such as two-party and multiparty 
teleportation H26] [J69l [J28l H70j HTB , and quantum key distribution fljS [139] EH ; they 
have been proposed for achieving one-way quantum computation with CV generalizations 
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of cluster states IU52L and in the multiparty setting they have been proven useful to solve 
Byzantine agreement II 1 731 . In this respect, one theoretical result of direct interest for 
the characterization of entanglement in Gaussian states is the qualitative and quantitative 
equivalence 01251 between the presence of bipartite (multipartite) entanglement in two-mode 
(TV-mode) fully symmetric Gaussian states shared as resources for a two-party teleportation 
experiment II 1261 11691 (TV-party teleportation network II 1281 11701 ). and the maximal fidelity 
of the protocol 111291 11301 . optimized over local single-mode unitary operations performed 
on the shared resource. In the special case of three-mode, pure GHZ/W states, this optimal 
fidelity is a monotonically increasing function of the residual contangle, providing the latter 
with a strong operational interpretation. Based on this equivalence, one can experimentally 
verify the promiscuous sharing structure of tripartite Gaussian entanglement in such states in 
terms of the success of two-party and three-party teleportation experiments [ 143]. 

Gaussian states are currently considered key resources to realize light-matter interfaced 
quantum communication networks. It has been experimentally demonstrated how a coherent 
state of light can be stored onto an atomic memory [174], or teleported to a distant atomic 
ensemble via a hybrid light-matter two-mode entangled Gaussian resource 11751 . 

Gaussian states play a prominent role in many-body physics, being ground and thermal 
states of harmonic lattice Hamiltonians [44 1. Entanglement entropy scaling in these systems 
has been shown to follow an area law 1 1761 11771 . In this context, entanglement distribution 
can be understood by resorting to the "matrix-product" framework, which results in an 
insightful characterization of the long-range correlation properties of some harmonic models 
[ 140 1 . Thermodynamical concepts have also been applied to the characterization of Gaussian 
entanglement: recently, a "microcanonical" measure over the second moments of pure 
Gaussian states under an energy constraint has been introduced |178|, and employed to 
investigate the statistical properties of the bipartite entanglement in such states. Under 
that measure, the distribution of entanglement concentrates around a finite value at the 
thermodynamical limit and, in general, the typical entanglement of Gaussian states with 
maximal energy E is not close to the maximum allowed by E. 

A rather recent field of research concerns the investigation of Gaussian states in a 
relativistic setting. The entanglement between the modes of a free scalar field from the 
perspective of observers in relative acceleration has been studied HI 791 11801 . The loss of 
entanglement due to the Unruh effect has been reinterpreted in the light of a redistribution of 
entanglement between accessible and unaccessible causally disconnected modes [ 1 80 1 . Such 
studies are of relevance in the context of the information loss paradox in black holes 111811 . in 
the general framework of relativistic quantum information 111 821 . 

From a fundamental point of view, some important conceptual and foundational 
problems include the operative interpretation and the possible applications of the relative 
entropy of entanglement for quantum states in infinite-dimensional Hilbert spaces 11831 . 
proving whether the Gaussian entanglement of formation coincides with the true entanglement 
of formation or not, elucidating the hierarchical structure of the distributed Gaussian 
entanglement between N modes, relative to all possible multipartitions k < N, and extending 
the quantitative investigation and exploitation of CV entanglement in a relativistic setting 
to the domain of squeezed and photon-augmented non Gaussian states 1 1 841 . Interest in 
non Gaussian entanglement is not restricted to foundational questions. In certain cases non 
Gaussian states may prove useful as resources in quantum information protocols, for instance 
for teleportation with non Gaussian mixed-state resources closely resembling Gaussian ones 
1 185|, quantum information networks with superpositions of odd photon number states 
[186|, and the experimental measurability of some quantifiers of entanglement such as the 
logarithmic negativity 118711 . It has been recently proven that there exist some particular 
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classes of non Gaussian squeezed states that can be produced and used as entangled resources 
within the standard Braunstein-Kimble teleportation protocol M126L and allow a sharp 
enhancement of the fidelity of teleportation, compared to Gaussian resources at fixed degree 
of squeezing 11 881 . 

We have seen how the monogamy constraint establishes a natural ordering and a 
hierarchy of entanglement of Gaussian states, that goes beyond the frustration effects that 
arise, for instance, in symmetric graphs 11411 . It is then natural to investigate how and in 
what form monogamy constraints arise in non Gaussian states of many-body CV systems, 
and, perhaps even more important, what is the structure of distributed entanglement for hybrid 
many -body systems composed of qubits, and/or qudits, interacting with single-mode or multi- 
mode fields 111891 . These investigations may be of special interest for the understanding 
of decoherence and entanglement degradation in different bath configurations, and the 
development of possible protection schemes 1 190|. From this point of view, some interesting 
hints come from a recent study of the ground-state entanglement in highly connected systems 
made of harmonic oscillators and spin- 1/2 systems [191|. We may expect that this area of 
research has more surprising, only apparently counterintuitive, results in store, besides the 
recent finding that, in analogy with finite-dimensional systems, independent oscillators can 
become entangled when coupled to a common environment [192] 

A fundamental achievement would of course be the complete understanding of 
entanglement sharing in a fully relativistic setting for general interacting systems with a 
general tensor product structure of individual Hilbert spaces of arbitrary dimension. In a 
non relativistic framework, the investigation of the structure of entanglement in hybrid CV- 
qubit systems is not only of conceptual importance, but it is relevant for applications as well. 
Here, we should at least mention a proposal for a quantum optical implementation of hybrid 
quantum computation, where qubit degrees of freedom for computation are combined with 
quantum continuous variables for communication 1119311 , and a suggested scheme of hybrid 
quantum repeaters for long-distance distribution of quantum entanglement based on dispersive 
interactions between coherent light with large average photon number and single, far-detuned 
atoms or semiconductor impurities in optical cavities [ 19411 . A hybrid C V memory realized by 
indirect interactions between different modes, mediated by qubits, has been recently shown 
to have very appealing features compared to pure-qubit quantum registers 11951 . 

It seems fitting to conclude this review by commenting on the intriguing possibility of 
observing CV entanglement at the interface between microscopic and macroscopic scales. 
In this context, it is encouraging that the existence of optomechanical entanglement between 
a macroscopic movable mirror and a cavity field has been theoretically demonstrated and 
predicted to be quite robust in realistic experimental situations, up to temperatures well in 
reach of current cryogenic technologies 111961 . This examples together with the few others 
mentioned above should suffice to convince the reader that continuous-variable entanglement, 
together with its applications in fundamental quantum mechanics and quantum information, is 
a very active and lively field of research, where more progress and new exciting developments 
may be expected in the near future. 
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